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Die Absorption in einem Elektronen-Gas-Gemisch 
Von Hermann Kober 
(Mit 1 Abbildung) 


Inhaltsiibersicht 


Der Einfluß der Stoßdämpfung auf die Polarisation eines unter der Wirkung 
eines elektrischen Wechselfeldes stehenden Elektronen-Gas-Gemisches wird auf 


' gastheoretischer Grundlage behandelt und der nichtlineare Effekt der Absorption 


in erster Näherung berechnet. 


Einleitung 


Für freie Elektronen, die sich in einem Gas befinden, gibt die Dispersions- 
theorie den Wert der Polarisation 


Ne 
= 


ot wenn ©, = Eei®t, 
oder in reeller Schreibweise 


2 
B,=— BER... (C08 ot sin wt), wenn ©, =E cos wt, (1) 
bzw. 
2 2 
=— (cos ot— sin t). wenn (=) (la) 


Dabei sind e und m Ladung und Masse des Elektrons, N die Zahl der Elektronen 
pro Volumeinheit und S die Stoßzahl, d. h. die mittlere Zahl der Stöße eines Elek- 
trons mit Gasmolekülen in der Zeiteinheit. 

H. A. Lorentz!) hatte im Dämpfungsglied obiger Formeln den Wert 2 S 
an Stelle von S errechnet, nach Salpeter?), Lassen?) und Krauß?) ist der ein- 
fache Wert S richtig. 

Die Betrachtungsweise aller genannten Autoren ist gekennzeichnet durch 
folgende vereinfachende Annahmen: 

a) Die Stoßzahl S wird als eine zeitlich konstante Größe angesehen, die ihren 
Wert unabhängig von der Phase der periodischen Elektronenbewegung, die sich 
der Wärmebewegung überlagert, bewahrt. 

Dies Verfahren wird näherungsweise angängig sein, solange die durch das Feld 
den Elektronen erteilten zusätzlichen Geschwindigkeiten klein bleiben gegenüber 


1) H. A. Lorentz, The theory of electrons, 1909, S. 306. 

2) R. Salpeter, Physik. Z. 14, 201 (1913). 

3) H. Lassen, Hochfrequenztechn. u. Elektroak. 28, 139 (1926). 
4) F. Krauß, Dissertation Köln 1950. 
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der mittleren thermischen Geschwindigkeit der Elektronen, d. h. solange E unter 
einer gewissen Grenze bleibt. Bei großen Feldstärken aber, wie siewnter Umständen 
in der Ionosphäre auftreten können, wird eine merkliche Modulation der mittleren 
Geschwindigkeit durch das Feld eintreten und infolgedessen auch die Stoßzahl 
im Rhythmus des Feldes schwanken. So entsteht durch die Stoßdämpfung ein 
nichtlineares elektrisches Verhalten des Elektronengases, das in der Ionosphäre 
den bekannten Kreuzmodulationseffekt (Luxemburgeffekt) zufolge hat. 


b) Die mittlere Stoßzahl S wird als Stoßwahrscheinlichkeit nicht nur für 
Elektronen aus der Gesamtheit aller — aus allen Geschwindigkeitsbereichen, — 
angesehen, sondern auch als Stoßwahrscheinlichkeit für Elektronen aus gewissen 
Auswahlen aus dieser Gesamtheit. Die Berechnung des Einflusses der Stöße er- 
folgt nämlich durch Erfassung der Zahl der seit einem gewissen Zeitpunkt unge- 
stoßen verbliebenen Elektronen und die Abnahme dieser ungestoßenen Elektronen 
wird durch das Gesetz e”S! beschrieben. Dabei wird nicht berücksichtigt, daß die 
seit einer gewissen Zeit ungestoßenen Elektronen eine Auswahl aus der Gesamtheit 
aller darstellen, in welcher im Laufe der Zeit die langsameren Elektronen immer 
mehr vorherrschen werden, weil die schnelleren eher zum Stoß kommen. Für diese 
Auswahl kann aber S nicht mehr die Stoßwahrscheinlichkeit sein. 


Im folgenden wird eine Näherungsberechnung der Polarisation auf der Grund- 
lage der Gastheorie durchgeführt, welche die beiden Vereinfachungen a) und b) 
nicht macht und damit vor allem einen Zugang zu der Nichtlinearität des elektri- 
schen Verhaltens des Elektronengases öffnet. Freilich muß eine solche Rechnung 
wegen der nun zu berücksichtigenden Geschwindigkeitsabhängigkeit der Stoß- 
wahrscheinlichkeit gewisse Annahmen über den Stoßmechanismus machen. Es 
wird die einfachste Vorstellung des elastischen kräftefreien Stoßes zugrunde gelegt. 


1. Verteilungsfunktion und Polarisation 


In dem unter der Wirkung des Feldes ©, = E cos w t stehenden Elektronengas 
ist die Geschwindigkeitsverteilungsfunktion f als Funktion der drei kartesischen 
Geschwindigkeitskomponenten «, v und w und der Zeit anzusetzen: f = f (u,v,w,t). 
Für E — 0 besteht Kugelsymmetrie im Geschwindigkeitsraum. Für E + 0 muß, 
da der die Kugelsymmetrie störende Feldeinfluß auf die z-Achse beschränkt ist, 
Symmetrie der Verteilung um die w-Achse des Geschwindigkeitsraumes bestehen, 
d.h. es muß f = f (c,d,t) sein, wo c = Yu? + v? + w? die Absolutgeschwindigkeit 
und ® der Winkel zwischen dem Geschwindigkeitsvektor ¢ = (u,v,w) und der 
z-Richtung ist. 


Deshalb kann man f nach den Kugelfunktionen des cos ® entwickeln: 


j= F, (ct) P, (cos). (2) 


Als Verteilungsfunktion muß f die Eigenschaft [fdr = 1 haben, wenn dr ein 


Element des Geschwindigkeitsraums ist und über diesen ganzen Raum integriert 


wird. Mit Ansatz (2) bedeutet das 


co xz 
2a f f SF, P, (cos #) ce de sin 3 di — 1. 
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Da 
wird 
J F, (c, t) ce? de = (3) 
Fiir die Polarisation gilt : 7 = N ew,, wobei über alle Elektronen in der Volum- 


einheit zu summieren ist, oder 
Nef fwd. (4) 


Da w = c cos v, erhält man mit (2) 


_ 2nNe> Sr. (c, node fr, (cos 9) sin cos dd. 


Es ist 
3fürn=1 


n +1 
J Pa (cos) sind cos P,, (x) = 0 fir n= 0,2,3... 


Also wird 
= J F, (ce, t) ede. (5) 


Für die Polarisation ist also nur der erste Koeffizient F, der Entwicklung von f 
maßgebend. 


2. Die Boltzu.annsche Gleichung und die nullte Näherungslösung 


Für die Zeitabhängigkeit der Verteilungsfunktion muß die Boltzmannsche 
Gleichung?) gelten. Sie läßt sich in unserm Fall schreiben als 


of ss Of eH 


da die Kraft auf die Masseneinheit des Elektrons = cos wt beträgt und z-Richtung 


hat. Die Gleichung stellt die zeitliche Änderung der Verteilungsdichte an einer 
Stelle des Geschwindigkeitsraumes durch zwei Glieder dar. Das erste Glied enthält 
die Dichteänderung durch einfaches Nachrücken von Elektronen an die betrach- 
tete Stelle des Geschwindigkeitsraumes infolge der durch das Feld verursachten 
Beschleunigung, der zweite Summand, das mit J abgekürzte Stoßintegral, enthält 
die Differenz des durch die Stöße verursachten Zugangs und Abgangs von Elek- 
tronen aus bzw. zu anderen Geschwindigkeiten. 

Vernachlässigen wir zunächst den Einfluß der Stöße vollständig, so erhalten 
wir für die nullte Näherung / die Gleichung 

(0) 
5) z. B. J. H. Jeans, Dynamische Theorie der Gase, 1926, 265ff. 
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mit der Lösung 


jo — (5 =. 7) ls u | (7) 
In dieser Verteilung ist gegenüber der Ruheverteilung (E = 0) die Geschwindig- 
keitskomponente w um den periodisch wechselnden Betrag Fan sin @ t verschoben, 
d.h. die Verteilung beschreibt einen Zustand, in dem die Gesamtheit der thermisch 
bewegten Elektronen zusätzlich mit der in 2- 
Richtung hin- und herschwingt. 


Mit (7) läßt sich aus (4) die Polarisation in nullter Näherung errechnen, sie er- 
gibt sich zu 
Ne&eE 


PB,=— mar cos wt, 


was mit (1) übereinstimmt, wenn S = 0 gesetzt wird. 


Mit Rücksicht auf späteren Gebrauch soll der erste Entwicklungskoeffizient 
F der nullten Näherungslösung / in seiner Abhängigkeit von c und ¢ dargestellt 
werden. 


Zunächst läßt sich in Polarkoordinaten und mit der Abkürzung a = sin wt 


schreiben als 


jo =( 2ac cos#) 
2a 


oder, entwickelt nach Potenzen von cos ?#, 


(Hat) 
(0) m 


Stellt man andererseits die Winkelabhängigkeit von / wieder durch Entwicklung 
nach Kugelfunktion dar: 


jo = 2 F® (c, t) P, (cos®), 
n= 


so ergibt sich auf Grund der Orthogonalitätsrelationen der Kugelfunktionen der 
erste Entwicklungskoeffizient zu 


+1 +1 
Fre — 35 1 P, (cos 9) d cos 8 = f® cos (cos d),. 


d.h. ausintegriert 
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3 
Zieht man auch noch den Faktor ¢ ##?" in die Entwicklung nach Potenzen von a, 


so erhalt man 


mit der Abkürzung 


(8) 


ez . 
a = ——- sinwt. 
mw 


3. Das Stoßintegral unter Annahme elastischer kräftefreier Stöße 


Um zu höheren Näherungen zu kommen, ist die Behandlung des Stoßintegrals J 
in der Boltzmannschen Gleichung nötig. Man erhält eine einfache Darstellung 
dieses Integrals, wenn man die Annahme des elastischen kräftefreien Stoßes ein- 
führt und die thermische Bewegung der Gasmoleküle vernachlässigt. 

Zunächst darf man von Stößen von Elektronen untereinander absehen, da sie 
wegen der kleinen Wirkungsquerschnitte — und in der Ionosphäre wegen der gegen 
die Gasdichte geringen Elektronendichte — gegenüber den Stößen von Elektronen 
mit Gasmolekülen, den gemischten Stößen, keine Rolle spielen. 

Für die gemischten Stöße nehmen wir an, daß die Stoßpartner sich wie ela- 
stische Kugeln verhalten und können dann bekanntlich das Stoßintegral schreiben 


J=o®N, [SSIS 9) ¢, cos O sin O dO de du, dv, dw. 
Hierin bedeuten 


o die Summe der Radien der beiden Stoßpartner, 

N, die Zahl der Gasmoleküle in der Volumeinheit, 

9: =9g(u,v,w,) die Geschwindigkeitsverteilungsfunktion der Gasmolekiile, 

c, der Betrag der Relativgeschwindigkeit der beiden Stoßpartner, 

© den spitzen Winkel zwischen der Zentrilinie des Stoßes und der Relativge- 
geschwindigkeit, 

e das Azimut der Zentrilinie in einem EEE dessen Achse 
mit der Relativgeschwindigkeit zusammenfällt. 

f = f(wW,v,w) bzw. gi = (uj, vj, wi) sind die Verteilungsfunktionen der 
Elektronen bzw. Moleküle an den Stellen u’, v', w' bzw. uj, vi, wy des Geschwin- 
digkeitsraums, an denen das Elektron bzw. sein Partner vor dem Stoß sein müssen, 
um durch den Stoß an die Stellen wu, v. w bzw. u,, v,, w, zu gelangen. 


Infolge des groBen Massenunterschiedes zwischen Molekiilen und Elektronen 
sind die thermischen Geschwindigkeiten der schweren Moleküle gering gegenüber 
den Elektronengeschwindigkeiten. Näherungsweise kann man deshalb die Tem- 
peraturbewegung des Gases vernachlässigen, d. h. die Moleküle als feste Hindernisse 
betrachten, an welchen die Elektronen gestoßen werden. 

Da dann g (u, v, w,) nur für u, = v, = w, = O nicht verschwindet und an dieser 
Stelle den Wert 1 hat und da ferner c, = c wird, erhält man 


aj/2 2% 


J=oN, ccos @ sin@ dO de. 
=0 ı= 
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Schreibt man in Polarkoordinaten f = f(c, ?, t), so wird f' = f (c, 8, t), weil beim 
elastischen Stoß am festen. Hindernis die Absolutgeschwindigkeit unverändert 


bleibt. Durch Entwicklung nach. 


Kugelfunktionen wird 


1-20 


fm =F. (c, t) P, (cos 8’). 


Die Abbildung zeigt die 
Beziehung zwischen ®’ und den 
3 übrigen Winkeln. Liegen der 

Vektor ¢ = (u,v,w) und die 
Zentrilinie 3 in der Zeichen- 
ebene, so liegt nach den Ge- 
setzen des elastischen Stoßes 
auch ¢’ = (u’,v’,w’) in dieser 
B+T-E Ebene und schließt mit c den 
Abb.1. Geschwindigkeiten vor und nach elasti- Winkel — 26 ein. 

schem Stoß am festen Hindernis Für den Winkel 9° zwischen 


ce’ und der nicht in der Zeichen- 
ebene liegenden z-Achse liefert der Cosinussatz 


cos = cos d cos (a— 20) + sind sin — 20) cos (®+n—e), 
wobei ® das Azimut der z-Achse in einem Polarkoordinatensystem mit c als Achse 
ist. 
Nach dem Additionstheorem der Kugelfunktionen ®) wird dann 


2 


f= = PR (cos 8) - P™ (— cos 20) eim(®+2- e), 


wobei pr die normierten zugeordneten Kugelfunktionen erster Art sind. 


Setzt man dies f’ in das Stoßintegral ein und führt zunächst die Integration 
über ¢ von 0 bis 277 durch, so verschwindet infolge der Periodizität der e-Funktion 
das Integral für alle Summanden mit m =+ 0 und für m = 0 wird der Integrand 
unabhängig von e. So erhält man 


n/2 
J=2no®N,c Pi F,, P,, (cos 8) J {P,, (— cos 2 @) — 1} cos O sin 0 dO. 


Das letzte Integral wird gleich 


fürn =0 


+1 


Damit erhalt man endgiiltig fiir das StoBintegral 


J=—no!N, P, (cos 9). (9) 


"s) Jahnke-Emde, Tafeln höherer Funktionen, 1948, S. 112. 
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-J 


4. Die erste Näherungslösung 


Zur weiteren Behandlung der Boltzmannschen' Gleichung entwickeln wir 
auch das zweite Glied 
öfeE 


— coswt 


nach Kugelfunktionen. Diese Entwicklung sei hier nicht ausgeführt (siehe auch 
die folgende Bemerkung), sie findet sich bei R. Lichtenstein’), der ähnliche 
Berechnungen für statische elektrische Felder durchgeführt hat. Sie ergibt — 


Ki n n(n —1) n+1 
1 
wobei F} = bedeutet. 


Mit dieser Entwicklung und dem Stoßintegral (9) läßt sich die Boltz mannsche 
Gleichung schreiben 


__ n(n—1) 1 
at Pp =— wt Fi P, tack 


Infolge der Orthogonalitätseigenschaften der Kugelfunktionen zerfällt diese 
Gleichung in das System 


— cos wt (Fi + 

usw. 


Um eine erste Näherungslösung zu gewinnen, benötigen wir nur die zweite Gleichung 
des Systems 10. Wir multiplizieren sie mit c® und integrieren über c von 0 bis oo. 
Dabei wird 
3 


oo oo co 
f Fed=F,e | —3 f Fuetde =——, 
6 0 0 4a 


weil c? F, auch für e =oo verschwinden muß und das letzte Integral nach (3) 
den Wert > hat. 
4x 


Ferner verschwinden die Glieder mit F,, weil 


=f f Fred =0. 
0 0 0 0 0 


7) R. Lichtenstein, Physik. Z. 39, 646 (1938). 


eim 
lert 
ach. 
) 
. 
den 
der 
die 
en- = 
Ge- 
ser 
den 
hen 
en- 
| 


8 Annalen der Physik. 6. Folge. Band 11. 1952 


W ir erhalten also 
ofr, ede =p eos fr 4.de. 


(Bemerkung: Gl. (11), die wir allein zu einer Näherungslösung benutzen werden, er- 


gibt sich auch ohne die Entwicklung von a nach Kugelfunktionen, wenn man die 


Boltzmannsche Gl. (6) mit w multipliziert und über den ganzen Geschwindigkeitsraum 
integriert: 


eE [ of 
feat cos wt [twat + wat 


co 
Darin ist f fwdt = An f Fc? de nach (4) und (5), ferner 
0 


u=-00 v=-0O 


und endlich mit (9) 
oo n/2 
ctde [ P, (cos) cos sind dd 
2=1 0 


[Fe de. 


4 20 eE 4n? 
Also wird Fre de = cos wi — o® N, [ Fıctde. 
Die Entwicklung von ai nach Kugelfunktionen war also für unsere Zwecke nicht 


nötig. Indessen kann das Gleichungssystem (10) eine Grundlage für höhere, hier nicht in 
Betracht gezogene Näherungen darstellen und sollte daher mitgeteilt werden.) 


co 
Gl. (11) legt folgende Näherungslösung für f F,c®denahe: Denkt man sich F, 
6 


nach o? entwickelt und vernachlässigt man Glieder mit o* und höheren Potenzen, 
so ist 


mit F® nach Gl. (8), welches wir aus der nullten Näherung berechnet hatten, 
Dann gilt in erster Näherung 


ar 3 eE 
J Pe de tn m no N J ctde (11a) 
Die Näherung besteht also darin, daß im Stoßglied die Verteilungsfunktion 
in nullter Näherung, d.h. so, als würden Stöße überhaupt nicht stattfinden, ein- 


gesetzt wird. Aus (8) erhält man durch Multiplikation mit c* und gliedweises 
Ausintegrieren 


“sich 
aus 


mit 


zu 


Dies 
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Dies in (lla) eingesetzt, ergibt mit (5) fiir die Polarisation, wenn wieder 


a= *—sinwt gesetzt wird 
mo 


70 3 


dt k T\mo 


Hier erscheint die Polarisation mit nichtlinearen Dämpfungsgliedern. In der 
Potenzreihe der sin wt haben zwei aufeinanderfolgende Glieder — abgesehen von 
dem immer kleiner werdenden Zahlenfaktor — das Koeffizientenverhältnis 
2 

(4) . Da die mittlere thermische Geschwindigkeit ¢ = 2 V — beträgt 
und die Amplitude der durch das Feld zusätzlich erzeugten Geschwindigkeit, 
wenn keine Stöße berücksichtigt werden, cp = = ist, wird das Koeffizienten- 


2 
verhältnis gleich S(% , also klein, solange die durch das Feld erzeugte Ge- 


schwindigkeit klein gegen die thermische Geschwindigkeit bleibt. 


In der Ionosphäre wird deshalb die Nichtlinearität vor allem bei Langwellen, 
bei denen cy groß wird, eine Rolle spielen können. Wählen wir w = 1,5. 10% 
(Langwelle), 7’ = 400 und setzen wir für die Feldstärke zunächst 10-® (das ist die 
Feldstärke in 100 km Entfernung von einem 150 kW-Sender bei kugelsymmetri- 
scher Abstrahlung in den leeren Raum), so ergibt sich mit k = 1,4 - 1071, m = 
9.102 und e = 5 - 10-10 für das genannte Verhältnis nur 2 - 10-3. Bedenkt man 
aber, daß in der Ionosphäre Feldstärken auftreten können, die um eine Zehner- 
potenz und mehr größer sind als es dem einfachen 1/r-Gesetz entspricht und daß 
die Feldstärke quadratisch in das Koeffizientenverhältnis eingeht, so wird es ver- 
ständlich, daß das sin? w t-Glied eine merkliche Größe erreichen kann. 


Berücksichtigt man nur die ersten beiden Glieder der sin w t-Reihe, so erhält 
man nach zweimaliger Integration aus (12) die Polarisation 


NeE akT 1 3 m /eE\*\ . 


@ 


— =) sin swt}. 


Die mittlere StoBzahl S bei ungestörter (E = 0) Elektronenbewegung ergibt 


‘sich unter unseren Voraussetzungen elastischer Stöße mit ruhenden Gasmolekülen 


aus 
S=N,o: [ fccos@ sin @ dO dedr 
me! 
mit f= (sen) e 
zu S=4N,0o® 
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Damit läßt sich die Polarisation schreiben 
NeE 48 3 m [eE\?\ . 
P,=— [(1 + ( ))sinot 


40kT\mo 
360 kT ot] 
2 
Für kleine Feldstärkewerte können die Glieder mit 7 2) vernachlässigt werden’ 


und es wird 


Der Vergleich mit (la) zeigt, daß im Dämpfungsglied statt S der Faktor 4/3 8 
auftritt. Wie in der Einleitung unter b) ausgeführt, ist der einfache Faktor S im 
Dämpfungsglied mit dem unerlaubten Gebrauch von S als Stoßwahrscheinlichkeit 
auch für gewisse Auswahlen aus der Gesamtheit aller Elektronen verknüpft, 
wird also nicht exakt richtig sein. Die Gültigkeit des Faktors 4/3 S nach unserer 
Rechnung ist andererseits durch die idealisierenden Voraussetzungen des elasti- 
schen Stoßes eingeschränkt. 


5. Die Polarisation bei großen Stoßzahlen 


2 
Wird die Stoßzahl so groß, daß (2) nicht mehr gegen 1 vernachlässigt werden 
kann, so enthält nach (1) die Polarisation den mit wachsendem S abnehmenden 
Faktor 


Eine entsprechende Abnahme der Polarisation mit wachsender Stoßzahl ist 
in unsern Resultaten (14) und (14a) nicht enthalten. Das liegt daran, daß wir 
im Näherungsansatz (lla), indem wir im Stoßglied die Verteilung f" bzw. F 
benutzten, uns auf in S lineare Glieder beschränkt hatten. 

Eine höhere Näherung, die auch die quadratischen Glieder der Stoßzahl S 
enthält, gewinnen wir, wenn wir im Stoßglied der Gl. (11) statt der Verteilung / 
nach (7) eine Verteilung 


(u? +0? + (w—b)*) 
(m) (15) 
mit i 
(sin ot) (15a) 
3 
ansetzen. 


Diese Verteilung liefert in (4) eingesetzt a eb und damit genau die 


Polarisation (14a). Sie stellt einen Zustand dar, in dem die thermisch bewegten 
Elektronen zusätzlich mit der Geschwindigkeit b pendeln; b ist aber die zur Nä- 
herungslösung (14a) gehörige mittlere aus dem Feld stammende Geschwindig- 
keit, die also durch Ansatz (15) den Elektronen jedes beliebigen Geschwindigkeits- 
bereichs zusätzlich zur thermischen Bewegung zugeschrieben wird. 

Geht man mit dieser Verteilung (15) in das Stoßglied der Gl. (11) ein, so unter- 
scheidet sich die Rechnung von der früheren nur dadurch, daß an Stelle von 


= ss sin w t jetzt b nach (15a) zu setzen ist. 


die 


Mit 


Fü 


wa: 
unc 
Ge: 


14) 
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Vernachläseigt man Glieder Sy, yT und höhere Potenzen und 


die gemischten Glieder von (= ya m0) ab, so erhält man 
kT\m 


B. = — [eos we (1— (5))-3 + 
_ 


Mit der gleichen Vernachlässigung von (2) läßt sich dies schreiben: 


%,=— {cos t— (1 sin w t 
mw? (1 + (42))) 3o 40k T\mo 
3o 


(16) 
1 m/eE 
—560£ 
Für den Fall kleiner Feldstärken ergibt sich 
Fe 


_ was bis auf die Faktoren 4/3 mit (1) übereinstimmt. 


Es ist selbstverständlich möglich, höhere Näherungsglieder zu berücksichtigen 
und statt 6 aus (15a) die der Polarisation (14) (statt (14a)) entsprechende mittlere 
Geschwindigkeit im StoBglied der Gl. (11) zu benutzen. 

Frechen-Buschbell, Ulrichstr. 97. 


(Bei der Redaktion eingegangen am 25. April 1952.) 
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Die optischen Konstanten von Gold und Silber 
im Ultraviolett 


Von Hermann Kober 


Inhaltsiibersicht 
Die optischen Konstanten dicker, im Vakuum aufgedampfter Gold- und 


Silberspiegel wurden mit einer sehr empfindlichen Methode, welche die Reflexionen | 


an 4 Spiegeln ausnützt, zwischen 550 mu und 250 mu bzw. 330 mu bestimmt. 


Die kürzlich erschienene Veröffentlichung von B. Seraphin!) über die Ultra- 
violettabsorption kolloidaler Goldlösungen, in welcher von mir gemessene, bisher 
unveröffentlichte optische Konstanten des Goldes zur Diskussion der Mieschen 
Theorie der Absorption von Kolloiden im Ultraviolett benutzt werden, gibt mir 


Veranlassung, meine 1939 durchgeführten Messungen an Gold und Silber kurz ' 


mitzuteilen, da neuere Messungen nicht vorzuliegen scheinen. 


1. Die Meßmethode 


Es wurden der Einfallswinkel g und das Azimut y der wiederhergestellten 
linearen Schwingung bei einer gewissen Phasenverschiebung 6 bestimmt. Zur 
Erhöhung der Meßgenauigkeit wurde dabei das an vier gleichen Metallspiegeln 
unter gleichen Einfallswinkeln sukzessive reflektierte Licht untersucht. Die ver- 
wendete Vier-Spiegel-Apparatur ist von R. Kretzmann?) beschrieben. Sie 
gestattete, die vier Spiegel mit großer Genauigkeit derart zu bewegen, daß der 
mit konstanter Richtung aus dem Polarisator eintretende Strahl an jedem der 
Spiegel unter gleichem, gemeinsam einstellbaren Einfallswinkel reflektiert wurde 
und danach wieder seine ursprüngliche Einfallsrichtung hatte, so daß die Analy- 
siereinrichtung ohne schwenkbaren Arm fest aufgebaut werden konnte. 

Als Analysator wurde ein Doppelbildprisma (Wollastonprisma) aus Quarz 
verwendet. Die beiden von ihm getrennten Komponenten wurden photographiert 
und die Schwärzungen photometrisch verglichen. 

Während bei der Methode des Doppelbildprismas, wenn sie nach einmaliger 
Reflexion angewandt wird, die Gleichheit der beiden Komponenten des mit 
seinen Hauptachsen unter 45° zur Einfallsebene gedrehten Doppelbildprismas die 


Phasenverschiebung = anzeigt, also nur beim Haupteinfallswinkel eintritt, ergibt 
sich bei Untersuchung des viermal reflektierten Lichts diese Gleichheit bei den 


? und —- = , d.h. bei den Phasenverschie- 
bungen 6 = + und! des einzelnen somit also bei vier verschie- 


1) B. Seraphin, Ann. Physik 10, 1 (1952). 
2) R.Kretzmann, Ann. Physik 37, 303 (1940). 
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37 
8 
Das dazugehörige Azimut y wurde dann in bekannter Weise bestimmt, indem 
das Doppelbildprisma um 45° gedreht, also mit seinen Hauptachsen parallel und 
senkrecht zur Einfallsebene gestellt wurde und das Polarisatorazimut so lange 
variiert wurde, bis wieder Gleichheit der beiden photographisch verglichenen 
Komponenten eintrat. Ist y das so gefundene Polarisatorazimut, so ergibt sich y 


aus tg p= Ves 


denen Einfallswinkeln. Die Messungen wurden stets bei 6 = durchgefiihrt. 


Als Lichtquelle wurde eine iiberlastete Kohlenbogenlampe benutzt, deren 
Licht durch einen UV-Monochrometer mit einem Abbeschen Prisma aus Quarz 
zerlegt wurde. Der Monochromatorspalt wurde stets so eingestellt, daß der Wellen- 
längenbereich des austretenden Lichts etwa 5 mu betrug. 

Die Empfindlichkeit der Messung war infolge der Benutzung der vierfachen 
Reflexion außerordentlich hoch. Der Fehler der Messung selbst betrug höchstens 
—+ 2 Minuten. Rechnet man dazu die Fehler: der sehr sorgfältig justierten Vier- 
spiegelapparatur mit + 4 Minuten, so kommt man auf höchstens + 6 Minuten 
Fehler in und y. Diese Fehler sind der Berechnung der Fehler in n und x in 
den letzten Spalten der Meßergebnistabellen zugrunde gelegt. Wie man sieht, 
gehen sie bei Silber viel stärker in die optischen Konstanten ein als bei Gold. 


. — Fehler, die in der Herstellung der Spiegel ihren Grund. haben, sind in diesen 


Angaben nicht berücksichtigt. 


2. Die Herstellung der Spiegel 


Die Spiegel wurden in einem Vakuum von ca. 6 - 10 Torr aufgedampft. Die 
Aufdampfanlage ist von R. Kretzmann (a.a. 0.) beschrieben. Die Spiegel 
wurden beim Aufdampfen zunächst mit einer Glimmerscheibe gegen die Schälchen, 
aus denen aufgedampft wurde, abgedeckt, damit die beim ersten Schmelzen des 
Metalls und Glühen der Schälchen etwa noch weggehenden Verunreinigungen 
sich nicht auf die Spiegel niederschlagen konnten. Nach dem ersten Glühen 
wurde dann die Glimmerscheibe umgeklappt (Bedienung von außen mittels 
Magneten) und die Spiegel so frei gegeben. Ferner wurde dafür gesorgt, daß sich 
stets Metall im Überfluß in den Schälchen befand, so daß nie bis zum letzten 
Rest verdampft werden mußte. Endlich wurde stets in Etappen aufgedampft, 
um eine zu große Erhitzung der Spiegel zu vermeiden. Die Dicke der Spiegel 
betrug bei Gold etwa 350 mu, bei Silber etwa 180 mu. 

Die vier zugleich aufgedampften Spiegel, die zu einer Messung verwendet 
wurden, waren im Sichtbaren innerhalb der Fehlergrenze der Messung ‚mit einem 
Polarisationsspektrometer stets optisch gleich. Durch sorgfältiges Konstant- 
halten aller Versuchsbedingungen beim Aufdampfen konnte einigermaßen Repro- 
duzierbarkeit der aufgedampften Spiegel erreicht werden. 


3. Die Meßergebnisse 


Die Tabellen 1 und 2 geben die gemessenen Werte für Gold und Silber. 

Infolge des etwas langwierigen photographischen Meßverfahrens, bei dem zu 
jedem Meßpunkt 9 Aufnahmen gemacht wurden, um die Stelle der Intensitäts- 
gleichheit der beiden Komponenten mit hinreichender Genauigkeit zu finden, 
dauerte die vollständige Messung einer Spiegelserie fast zwei Tage. Zum Ende 
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Tabelle 1. Gold = gr 
[77 y u x b= du dx 
550 74°25’ | 39°05’ | 0,5653 3,668 2,070 0,012 0,064 
500 71°47’ | 34°50’ | 1,006 1,690 | 1,700 
450 72°31’ | 28°54’ | 1,347 1,201 | 1,618 
400 72°40° | 27°17’ | 1,487 1,073 1,596 0,014 0,007 
370 72°18’ | 27°63’ | 1,411 1,111 1,567 
340 72°15’ | 27°13’ | 1,458 1,058 | 1,543 
310 72°11° | 27°36’ | 1,421 1,089 | 1,548 0,016 0,008 
280 71°19’ | 28°15’ | 1,284 1,165 | 1,495 
250 69°30’ | 28°07’ | 1,133 1,201 | 1,361 | 0,018 0,009 

3 

Tabelle 2. Silber 
| y u k=ux du dx 
550 | 80°02’ | 44°05’ | 0,1345 | 26,85 3,612 0,015 2,92 
500 | 78°45’ | 44°05’ | 0,1194 | 26,33 3,144 
450 | 77°00° | 44°02’ | 0,1097 | 24,11 2,646 
400 | 74°16’ | 43°44’ | 0,1205 | 17,18 2,022 0,009 1,35 
370 71°42’ | 43°12’ | 0,1501 | 11,10 1,665 
340 66°45’ | 42°00’ | 0,2103 | 5,207 1,095 
330 | 62°55’ | 41°02’ | 0,2640 | 2,888 0,7624 | 0,004 0,074 


der Messung wurde stets die erste Aufnahme vom Vortag wiederholt. Sie ergab 
bei Goldmessungen auf die Minute das gleiche Resultat. 

Anders bei Silber. Hier sanken im Lauf von 2 Tagen Einfallswinkel und 
Azimut um 30—40 Minuten ab. Genauere Priifung dieser Erscheinung zeigte, 
daß dieses Absinken der g- und y-Werte in den ersten zwei Stunden nach Heraus- 
nehmen der Spiegel aus dem Vakuum nur etwa 6—8 Minuten betrug. Es wurden 
deshalb mit einer Spiegelserie jeweils nur etwa 2 Stunden lang beschleunigte 
Messungen angestellt und fiir weitere Messungen neu aufgedampfte Spiegel ver- 
wendet, wobei durch Kontrollmessungen an gemeinsamen Meßpunkten die Gleich- 
heit der Spiegel verschiedener Serien bis auf wenige Minuten sichergestellt wurde. 
So sind die Ergebnisse der Tabelle 2 aus verschiedenen Serien von Silberspiegeln 
zusammengestellt. 

Bei dem Effekt des Absinkens der g- und y-Werte handelt es sich offenbar 
nicht um ein Altern infolge innerer Strukturänderungen der Schicht, sondern um 
Einflüsse der Luft oder von Fremdgasen; denn eine Spiegelserie, die 3 Tage lang 
im Vakuum, in dem sie hergestellt worden war, aufbewahrt wurde, zeigte Ein- 
fallswinkel und Azimute, die bis auf 5 Minuten mit denen frisch aufgedampfter 
Spiegel übereinstimmten. 


Frechen-Buschbell, Ulrichstr. 97. 


(Bei der Redaktion eingegangen am 15. April 1952.) 
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Über den Einfluß der Korngröße auf die Koerzitivkraft 
Von A. Mager 


Inhaltsübersicht 


Mit einem stark vereinfachten Modell der B>zirksstruktur läßt sich der Ein- 
fluß der Korngröße auf die Koerzitivkraft nach Größe und Gang richtig angegeben. 


Von T. D. Yensen und N. A. Ziegler?) wurde in gründlichen experimentellen 
Arbeiten u. a. auch der Einfluß der Korngröße (Querschnitt Q) auf die Koerzitiv- 
kraft (H) untersucht. Spätere Untersuchungen von W. E. Ruder?) sowie von 
O. Dahl, F. Pawlek und J. Pfaffenberger!) bestätigen die von Yensen und 
Ziegler empirisch aufgestellte Beziehung für 7, (Korngrößenanteil): 

H, (Oe) = & — 
Eine einfache Deutung ergibt sich, wenn man die Kristallkörner als Träger von 
Ummagnetisierungskeimen ansieht, deren Durchmesser durch die Abmessung 
dieser Körner bestimmt wird. Es gilt dann für die Startfeldstärke H, des unbe- 
schränkten Längenwachstums dieser Keime nach einer Beziehung von W. Döring?) 


_3ay 
worin y die spezifische Wandenergie der Blochwand (180°-Wand) und J, die Sätti- 
gungsmagnetisierung ist. Diese Beziehung wird im wesentlichen auch für kompli- 
ziertere Blochwandstrukturen erhalten bleiben infolge der Kopplung der Wand- 
verschiebungen durch die Streufeldenergie. 
Wenn man über alle Raumrichtungen mittelt, so ergibt sich noch ein Faktor 

von etwa 3/2). 

Mit den Werten yısoore © 1,5 erg/cm? 5) und J, = 1700 EME bei Fe folgt damit 
fir H, 
3 32-1,5 3,1-10-8 


H,[0e] 4-1700d [cm] d [cm] 


in relativ guter Übereinstimmung mit der von Yensen und Ziegler aufgestellten 
Beziehung. 


1) Nach F. Pawlek, ,,Magnetische Werkstoffe‘ Springer-Verlag, Berlin-Göttingen- 
Heidelberg 1952. 

2) W. E. Ruder, Trans. Amer. Soc. Metals 22, 1120/41 (1934). 

3) W. Döring, Z. Physik 108, 137 (1938). 

4) R. Becker u. W. Döring, ,,Ferromagnetismus“, Springer Verlag, Berlin 1939. 

5) Nach M. Kersten, Z. Physik 124, 714—741 (1947) aus 2 yo0°. 
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Ein gewisser Einfluß der Blechdicke auf die Koerzitivkraft kann mit dem glei- 
chen Effekt erklärt werden. Die Messungen von O. Dahl, F.Pawleck und J, 
Pfaffenberger!) zeigen jedoch besonders für sehr kleine Blechdicken eine wesent- 
lich stärkere Abhängigkeit, was nur durch ein grundsätzlich anderes Verhalten 
(höhere Koerzitivkraft) der Randzonen der Bleche nach F. Feldtkeller®) erklärt 
werden kann. 

Als weiterer Hinweis auf diesen Sachverhalt muß die Tatsache gedeutet werden, 
daß die von Sixtus und Tonks’) untersuchten Ummagnetisierungskeime sich 
meist im Inneren des Drahtes ausbilden, obwohl die Lage an der Oberfläche energe- 
tisch günstiger wäre wegen der Einsparung von Wandenergie. 


6) F. Feldtkeller, Fernmeldetechn. Z. 2, 9—14 (1949). 
*) K. J. Sixtus u. L. Tonks, Physic. Rev. 87, 93a (1931) u. weitere Arbeiten. 


Jena, Forschungsinstitut für magnetische Werkstoffe. 


(Bei der Redaktion eingegangen am 19. Juli 1952.) 
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Innere Feldemission aus lokalisierten Fehlstellen 
isolierender Kristalle 


Von Walter Franz 
(Mit 4 Abbildungen) 


Inhaltsübersicht 


Ein in einer Fehlstelle eingefangenes Elektron kann durch ein äußeres elek- 
trisches Feld mittels wellenmechanischen Tunneleffekts (innerer Feldemission) 
befreit werden. Um die zeitliche Wahrscheinlichkeit dieses Vorganges zu berechnen, 
wird die Fehlstelle als kugelförmiger Potentialtopf idealisiert, von der Anisotropie 
des Kristalls abgesehen und überdies die Fehlstelle als flach angenommen. Ver- 
nachlässigt man die (praktisch bedeutungslosen) Potentialdifferenzen, welche das 
Feld im Inneren und in der unmittelbaren Umgebung des Topfes hervorruft, so 
kann man die Wellenfunktion des Elektrons ohne weitere Vernachlässigungen 
durch ein handliches Integral darstellen, aus welchem sich die Austrittswahr- 
scheinlichkeit errechnen läßt. Das Ergebnis zeigt, daß die Emission aus Fehl- 
stellen nur dann merklich sein kann, wenn die Bindungsenergie unter etwa !/, Volt 
liegt; im Gegensatz dazu ist Feldemission aus dem Valenzband bei Energielücken 
bis 1 Volt möglich. 


Einleitung. Berechnung der Dämpfungsfaktoren 


Von Feldemission spricht man, wenn ein elektrisches Feld den Austritt von 
Elektronen mittels wellenmechanischen Tunneleffekts ermöglicht. Für den klassi- 
schen Fall der Emission aus einer kalten Metalloberfläche zeichnen wir in Abb. 1 
den Potentialverlauf ohne (gestrichelt) 
und mit Feld (ausgezogen). Die IE 
Metallelektronen liegen bei der abso- . 
luten Temperatur 7’ = 0 energetisch , A 
alle unterhalb der Fer mischen Grenz- 
energie & und werden durch die Tr 
Potentialschwelle von der Höhe &+ A 
(A = Austrittsarbeit) am Verlassen Abb. 1. Äußere Feldethission aus einer 
der Oberfläche gehindert, da der Metalloberfläche ins Vakuum 
Energie £ außerhalb des Metalls eine 
imaginäre Wellenzahl entspricht. Durch ein äußeres Feld F wird im Abstande 
A/eF von der Oberfläche ein Gebiet geschaffen, in welchem zu E = £ wieder 
reelle Wellenzahlen gehören, so daß dort eine auslaufende Elektronenwelle mög- 
lich ist. Das Zwischengebiet imaginärer Wellenzahlen wird durch Tunneleffekt 
überwunden, wobei aber die Welle stark gedämpft wird, und zwar ist der 
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Dampfungsfaktor der Intensität 


AleF 
D=exp(—21m f (1) 
6 
Fiir die Wellenzahl gilt die Beziehung 
h? k? = 2m (E— V) = —2m(A—e Fz). (2) 


Darin bedeutet (= £) die Gesamtenergie und V die potentielle Energie des be- 
trachteten Elektrons. Mit (2) ergibt sich aus (1) die bekannte Formel!) 


D,= exp(— 4) (3) 


Die Feldemission aus der kalten Oberfläche wollen wir als äußere Feldemission 


bezeichnen, im Gegensatz zu den verschiedenen Fällen innerer Feldemission, 


k 2 
Valenzband) Leitungsband 
verborenes Gebier 


Abb. 2. Innere Feldemission aus dem Valenz- Abb. 3. Wellenzahl als Funktion der 
band kinetischen Energie 


welche sich im Innern von Isolatoren abspielen können (da hierbei keine Bild- 


kraft auftritt, haben wir auch bei dem einleitenden Beispiel der äußeren Feld- 
emission von ihr abgesehen). 

Im Isolator haben wir ein vollbesetztes Valenzband mit reellen Wellenzahlen, 
welches durch ein ,,verbotenes“ Energiegebiet imaginärer Wellenzahlen von dem 
Leitungsband getrennt ist, in dem die Wellenzahlen wieder reellsind. Die Energie- 
lücke habe die Höhe J ( = Ionisierungsenergie). Durch ein äußeres Feld werden 
die Energiebänder gekippt (s. Abb.2), so daß ein Elektron des Valenzbandes nach 
Durchlaufen eines Tunnels von der Länge J/e F ins Leitungsband eintreten kann. 
Die Durchlässigkeit errechnet sich analog zu (1), doch ist der Zusammenhang 
zwischen Energie und Koordinate komplizierter als (2); am unteren Ende des 
Leitungsbandes ist kinetische Energie wie Wellenzahl (nach geeigneter Reduktion) 
gleich null, in der Energielücke ist k? negativ, wird jedoch am oberen Rande des 
Valenzbandes wieder null und dann positiv. Wir setzen einen parabolischen Ver- 
lauf für die Abhängigkeit zwischen kinetischer Energie 7 — £&—V und dem 
Quadrat k? der Wellenzahl an, wie er aus jedem Modell im Limes schwacher 
Störung durch das Gitter (d.h. kleiner Energielücke) folgt (s. Abb. 3): 


2k? = 2m* T (1 + T/I). (4) 


Darin ist die multiplikative Konstante m* die effektive Masse am unteren 
Rande des Leitungsbandes. Für ein Elektron, welches gerade soviel Energie be- 


1) Siehe z. B. H. Fröhlich, Elektronentheorie der Metalle, Berlin 1936, S. 88, Gl. (9). 
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sitzt, daß es bei x = 0 ins Leitungsband eintreten kann, ist 7 —eFx. Die Durch- 
lässigkeit errechnet sich zu 


0 

a V2m* 

= I = 

D, exp ( 21m Kae) (5) 

In dem betrachteten Fall kleiner Energieliicken hat man zwischen m*, J und 

Gitterkonstante a den Zusammenhang (s. H. Fröhlich, l.c. 8.45, Gl. (24)) 
2m? a? I 


(6) 
womit sich aus (5) die Zenersche*) Formel ergibt: 
m a I? 

D, = exp (- incr) . (7) 

Neben der inneren Feldemission aus dem Valenzband kommt auch 

die Emission aus einem im verbotenen Energiegebiet gelegenen Störniveau 

in Betracht, welches um den Betrag AE unter dem Rande des Leitungsbandes 

gelegen sei. Zur Berechnung des Dämpfungsfaktors haben wir auch hier Formel (5) 

zu verwenden, jedoch die untere Grenze der Integration nach —AZ/eF zu ver- 
legen. Das Ergebnis ist nach elementarer Integration 

* 
Ds = exp [- {arcos (1— )- 2(1— . (8) 


Für flach liegende Störniveaus (AH < 1/2) geht dies über in 


4 V2n m* 3 1 3 

DE AB "). (8a) 
Man ersieht aus dem Vergleich von (8) bzw. (8a) mit (5), (7), daß man erhebliche 
Fehler begeht, wenn man die Zenersche Formel auf die Emission aus Störstellen 
anwendet, indem man einfach J durch AZ ersetzt; dies kann nicht verwundern, 
da ja die Emission aus dem Valenzband dadurch erleichtert wird, daß k am oberen 
Rande dieses Bandes dem Wert null zustrebt, was auf dem Niveau einer Stör- 
stelle nicht der Fall ist: 

Wegen seiner rapiden Veränderlichkeit ist Fe exponentielle Dämpfungs- 
faktor allein für die Größenordnung der Feldstärken maßgebend, bei welchen 
Emission einsetzt. Für quantitative Betrachtungen benötigt man jedoch den 
vollen Ausdruck der zeitlichen Emissionswahrscheinlichkeit. Für die Emission 
aus dem Valenzband ergibt sich dieser einfach aus der Frequenz, mit welcher ein 
Valenzelektron zwischen den Rändern des Valenzbandes hin- und herpendelt, 
zu (s. H. Fröhlich, ].c. S. 47) 


Ds exp(— 


Ww, = (9) 


Die Berechnung der Emissionswahrscheinlichkeit aus einer lokalisierten Fehl- 
stelle ist der Inhalt der folgenden Rechnung. 


$1. Darstellung einer lokalisierten Fehlstelle als Potentialtopf 
Spaltet man von der Wellenfunktion eines Elektrons im. idealen Kristall- 
gitter entsprechend dem Blochschen Ansatz eine im Gitter periodische Modu- 
2%, Cl. Zener, Proc. Roy. Soc. London (A) 145, 523 (1934). 
QF 
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lation ab, so verbleibt eine Lösung der Wellengleichung 
(A+k)p=0. (10) 
Die Wellenzahl & hängt darin, wenn wir das Gitter als isotrop behandeln, in der 
einleitend geschilderten allgemeineren Form von der Energie ab. Beschränken 
wir uns auf die Betrachtung des 
Energiegebietes dicht ober- und 
unterhalb des unteren Randes des 
Leitungsbandes, dann können wir 
näherungsweise setzen 


k? = 2m*(E—V). (11) 


Wir stellen jetzt eine flache Fehl- 
stelle im Kristall dar als einen 
kugelförmigen Potentialtopf' vom 
Radius R, in dessen Inneren das 
Potential um einen konstanten Be- 
Abb. 4. Potentialtopf und Ladungsverteilung trag gegenüber dem umgebenden 
für xR = 27/3 (ky R = 22/3) ungestörten Kristall erniedrigt sei 
(Abb. 4). Die Tiefe des Topfes und 
gleichzeitig die energetische Lage des darin gebundenen Elektrons charakteri- 
sieren wir durch die Angabe der Wellenzahl &, im Inneren und der komplexen 
Wellenzahl ix im Äußeren des Topfes. Die Wellenfunktion des Elektrons können 
wir dann ansetzen: 2 


<P 
or 


y- (12) 
fürr>R. 
xr 
x hängt mit der Bindungsenergie AZ zusammen vermittels (11): 
x? 2m* AL. (13) 


Der Koeffizient A ergibt sich aus der Bedingung des stetigen Anschlusses von y 
und @y/ér bei r = R zu 
Ve+k 
Für die Wellenzahl k, des Grundzustandes erhält man gleichzeitig die transzen- 
dente Gleichung 


(14) 


tank, R = — kulx; 5 <kR<a. (15) 

Unsere Eigenfunktion ist nicht auf 1 normiert, vielmehr folgt aus (12), (13), (14) 
2n 

J (1 + R). (16) 


Das Elektron ist nur dann im Topf eingeschlossen, wenn «R>1; fiir xR <1 
dagegen hält es sich hauptsächlich in der Umgebung des Topfes auf. Einen mitt- 
leren Fall (xR = 22/3 - /3) zeigt Abb. 4. Wir werden später sehen, daß xR <1 
der Normalfall ist (s. Gl. (41)), somit das Elektron eines F-Zentrums sich nicht am 
Ort der Anionenlücke aufhält, sondern in der Umgebung vagabundiert, was der 
üblichen Auffassung durchaus entspricht. 
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Fiir die kommende Rechnung entwickeln wir die Wellenfunktion des AuBen- 
raums mittels einer Integraldarstellung von Sommerfeld?) nach Kegelwellen: 


oo 
=A f dy Sin x Jo(% o Sin x) Griz, (17) 
0 


+e (18) 


Darin ist 
der Abstand von der x-Achse. 


§ 2. Kegelwellen im homogenen elektrischen Feld 


Wir denken uns ein homogenes elektrisches Feld angelegt, welches so klein ist, 
daß sein Einfluß auf die Verhältnisse im und am Rande des Topfes vernachlässigt 
werden darf. Wir können dann an der Topfoberfläche die Integraldarstellung (17) 
verwenden; um auch entfernte Stellen zu erfassen, müssen wir an den Kegel- 
wellen diejenigen Modifikationen anbringen, welche das äußere Feld bewirkt. 
Die Wellengleichung, welche nach (10) und (11) jetzt lautet 


(4+ (19) 
läßt die gesuchten Kegelwellen sofort separieren in der Gestalt 
u = I (ko 0) (2). (20) 
Für den von x abhängigen Faktor erhält man die Gleichung 


Wollen wir unter % eine der Kegelwellen des Integranden von (17), modifiziert 
durch das Feld, verstehen, dann haben wir zu setzen: 


k, = x Gin x; v(0) = 1. (22) 


Wir führen in (21) als neue Variable ein 


2m* e F\'/s hi? (x? + 
(«— (23) 
und haben dann 
b+tv=0 (24) 
mit der Lösung (s. Magnus-Oberhettinger, l.c. 8. 26 u.) 
ts). (25) 


Unter den verschiedenen Zylinderfunktionen Zı,, haben wir die dem Problem an- 
gemessenen auszuwählen. Für sehr große ¢ ist die auslaufende Welle zu nehmen. 
also die Hankelfunktion erster Art. Allerdings können wir diese nur für große 
positive x unbesehen verwenden; der singuläre Punkt t = 0 der Zylinderfunk- 
tionen ist nämlich, wie (23) zeigt, im Gebiet <> 0 enthalten und erfordert eine 
gesonderte Betrachtung. Wir haben jedenfalls 


v= für (26) 


3) Siehe W. Magnus und F. Oberhettinger, Formeln und Sätze für die speziellen 
Funktionen der Mathematischen Physik, 2. Aufl., Berlin 1948, S. 48u. 
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Der Koeffizient B ist später zu bestimmen. — Für ¢ <0 können wir als vollstän- 
digen Satz von Zylinderfunktionen die beiden Hankelfunktionen vom Argument 
; %|t)/* heranziehen. Wir müssen sie so linear kombinieren, daß bei t = 0 die 
Stetigkeit von v und seiner Ableitung gewahrt wird. Durch Entwicklung um 
t = 0 erhält man leicht 


-in/6 1/, -i n/3° 
Somit haben wir zu setzen 


v= BY\t für (27) 


xz =. 0. 


Da diese Formel bei x = 0 gültig ist, gestattet sie mittels (22) die Bestimmung 
von B: 
B (28) 


Darin ist i, derjenige Wert von t, welcher sich für x — 0 ergibt: 


m* (x2 + ko) (29) 


2m*eF 
§ 3. Austritt aus dem Potentialtopf 


Indem wir die Kegelwellen (20) unter Benutzung von (26), (27), (28) in die 
Integraldarstellung (17) einfiihren, erhalten wir die folgende fiir > 0 im AuBeren 
des Topfes giiltige Darstellung fiir die Wellenfunktion des Fehlstellen-Elektrons: 


Z 
H? i |to 


an VL 

A Sin J, (xo Sin x) Va 
0 

mit 


2= ( für t > 0; — Hi für t<0; (31) 


* Is AE 
Bei nicht zu großen Feldern wird die iets dadurch schwach gehalten, 
daß t, von großem Betrage und daher die im Nenner von (30) enthaltene Hankel- 
Funktion sehr groß ist. Wir untersuchen nur diesen Fall schwacher Feldemission 
(bei starker Feldemission wird die Fehlstelle sofort entleert und damit eine quan- 
titative Durchrechnung überflüssig) und dürfen deshalb die Hankel-Funktion 
des Nenners durch ihre asymptotische Darstellung ersetzen. Betrachten wir 
weiterhin die Verhältnisse an einer Fläche x = const> A E/e F, so wird t> 1, 
abgesehen von extrem großen Werten von y, deren Beitrag wegen des exponen- 
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tiellen Anwachsens des Nenners sowieso vernachlässigbar ist. Wir dürfen deshalb 
für <> A E/e F Formel (30) ersetzen durch 


dz Gin Sin x) V Gof x 


er ( (2-42 cop (4 4E 


(33) 


Da der Exponent, wie vorausgesetzt, einen negativen Realteil von groBem Betrage 
besitzt (f,> 1!), liefern nur Gebiete kleiner x beträchtliche Beiträge zum Integral, 
und wir dürfen ohne merklichen Fehler Gin y durch x, Cojyx durch 1 + x?/2, 


außerhalb des Exponenten sogar durch 1 ersetzen. Dies ergibt: 


Darin bedeutet X das Weber-Soninesche Integral (s. Magnus-Oberhettinger, 
l.c. 8. 49): 


ı te 
(35) 


Als Gesamtausbeute der Feldemission berechnen wir die Strömung senkrecht 
zu der angenommenen weit entfernten Ebene x = const. Die Stromdichte ist 


h (ver dy V2 m* e Fx 


) dx m (36) 
me 
= m* ° 


Die zeitliche Wahrscheinlichkeit w, dafür, daß ein ee befreit 
wird, ist offenbar 


2x S jede 
Ws = . (37) 
Aus 66), (16), (14) erhalten wir schließlich 
1 
Für kleine xR (s. hierzu (41)) geht dies unter Benutzung von (6) über in 
aeF m ı/I heF 
Ws N oh AE Ds == 4maVIAB Dg. (38a) 


Dg ist genau der in der Einleitung berechnete Dämpfungsfaktor (8a). Wir haben 
dort gesehen, daß für tiefe Störstellen (AH nicht klein gegen J/2) statt dessen 
Formel (8) zu treten hat. Den Koeffizienten in (38) haben wir zwar unter der 
Voraussetzung flacher Störstellen berechnet, können ihn aber unbedenklich auch 
für tiefe Störstellen verwenden, da der dabei begangene Fehler neben der rapiden 
Feldstärkenabhängigkeit von D, und dem durch sie bedingten Fehler bei der 
Anwendung der Formel bedeutungslos bleibt. 
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Schlu8. Numerische Betrachtungen 


Um die praktischen Konsequenzen der abgeleiteten Formeln zu sehen, führen 
wir die Zahlenwerte der universellen Konstanten ein. J und AE messen wir in 
eVolt, R und a in Ä-Einheiten und F in Volt/cm und haben dann nach (9) und (7): 


w, = 2,42: 18a F exp(—3,28 - 10" 7) (39) 

Nach (38), (8a), (6), (13) ergibt sich: b 
F a 

exp ( 5,57: 100 ) (40) 
mit 

xR = V1 - AE; tank, R=— (41) 

i 


Man erkennt aus (41), daß xR nur bei besonders hohen Werten der Größen a, R, I, 
AE den Betrag 1 erreicht, während k, R ja größer als 2/2 ist. Das bedeutet, daß 
die ersten beiden Brüche in (40) sich gewöhnlich nur wenig von 1 unterscheiden. 
Aus den Exponenten von (39) und (40) ersieht man, daß w, und w, nur dann eine 
bemerkbare Größenordnung erreichen, wenn F ~ 10% (Volt/em) ist, womit die 
Koeffizienten der Exponentialfunktion in beiden Formeln rund 1012 bis 101% 
werden. w-Werte, von etwa 1/sec erhält: man daher, wenn die Exponenten gleich 
12 bis 13 mal In 10 sind, also rund 30. Emission aus Fehlstellen wird sich etwa bei 
Austrittswahrscheinlichkeiten dieser Größenordnung bemerkbar machen. Die 
Emission aus dem Valenzband dagegen kann wegen der um einen Faktor 105 bis 
10% größeren Dichte der Valenzelektronen schon bei einem um 10 Zehnerpotenzen 
kleineren w feststellbar sein. Daher erhalten wir für die Feldstärken F, bzw. F,, 
bei welchen Emission aus dem Valenzband bzw. den Störstellen merkbar wird, 
roh die folgenden Abschätzungen®) 


F, + 0,5-108aI?; Fs 2: 105a (AB), (42) 


Dies besagt, daß Feldemission aus dem Valenzband nur bei recht kleinen J in 
Frage kommt (nicht größer als etwa 1eVolt); andernfalls liegt F, über der Durch- 
schlagsfeldstärke. Die Feldemission aus Störstellen kann wegen des größeren 
Zahlenfaktors erst bei Bindungsenergien AZ unterhalb etwa 1/, Volt in Erschei- 
nung treten. 


4) Die Methoden einer genaueren Berechnung finden sich in W. Franz, Naturwiss. 
89, 18 (1952); Z. Physik 132, 285 (1952). 


Münster i. W., Institut für theoretische Physik der Universität. 


(Bei der Redaktion eingegangen am 5. Juni 1952.) 
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Strahlungskorrekturen zur Mesonerzeugung 
durch Photonen 


Von H. Lehmann 
(Mit 3 Abbildungen) 


Inhaltsiibersicht 


Es wird die zweite Näherung hinsichtlich der Meson-Nukleon-Kopplung be- 
rechnet fiir den Wirkungsquerschnitt der Erzeugung eines positiven Mesons durch 
Wechselwirkung eines Photons beliebiger Energie mit einem Proton. Zugrunde 
gelegt ist eine geladene pseudoskalare Mesontheorie mit pseudoskalarer Kopp- 
lung. Der Einfluß der Korrektur auf den differentiellen und totalen Wirkungs- 
querschnitt wird betrachtet. Im letzteren Fall ist (mit c= 0, + 0,) die Korrektur o, 


für alle Energien < 0,1 In hie 


I. Einleitung 


Obwohl auch heute noch zur feldtheoretischen Behandlung von Nukleon- 
Mesonprozessen im allgemeinen als Näherungsverfahren eine Entwicklung nach 
der Kopplungskonstanten benutzt wird, liegen nur wenige Rechnungen üder den 
Einfluß der zweiten nichtverschwindenden Näherung vor. Nur die renormali- 
sierbaren Theorien gestatten eine derartige Untersuchung; in erster Linie sind hier 
pseudoskalare Mesonen mit pseudoskalarer Kopplung von Interesse. Wir beziehen 
uns im folgenden auf diesen Fall. Betrachtet wurden unseres Wissens bisher die 
Korrekturen zur Streuung von Mesonen an Nukleonen!), sowie zur Erzeugung 
von Mesonen durch Photonen?); beides nur im Grenzfall kleiner Energien der 
beteiligten Teilchen. Es ist bekannt, daß die Zulässigkeit der erwähnten Näherungs- 
methode ernsten Zweifeln unterliegt; die Größe der Kopplungskonstanten läßt 
sich nicht mit einiger Genauigkeit angeben. Angesichts der Tatsache, daß bisher 
kein anderes konsequentes Verfahren entwickelt worden ist, scheint es dennoch 
wünschenswert, genauere Aussagen über den Beitrag der zweiten Näherung zu 
gewinnen. 

Wir betrachten deshalb fiir einen bestimmten Prozeß die Korrekturen zum 
Wirkungsquerschnitt. Die Rechnung wird für alle Energien durchgeführt. Die 
Photon-Mesonerzeugung stellt wohl den einzigen Prozeß dar, für den dies ohne 
größeren Aufwand möglich ist. 

Wir beziehen uns auf die Erzeugung eines 2+-Mesons in einer reinen geladenen 
Theorie. Wie erwähnt, sind die Korrekturen in der Nähe der Energieschwelle 
schon von Brückner?) untersucht worden; unsere Ergebnisse stimmen bei 
Spezialisierung hierauf jedoch nicht mit den dortigen Resultaten überein. Da in (2) 


1) Ashkin, Simon u. Marshak, Prog. Theor. Phys. 5, 634 (1950). 
2) K. A. Brückner, Physic. Rev. 79, 641 (1950). 
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keine näheren Angaben über den Gang der Rechnung gemacht werden, läßt sich 
die Ursache dieser Unterschiede nicht ohne weiteres angeben. 

Bei dem betrachteten Prozeß liegt eine Kopplung von drei Feldern vor. Die 
Renormalisierungstheorie hierfür ist u. a. von Ward?) behandelt worden. 

In der Bezeichnungsweise schließen wir uns den Arbeiten von Feynman‘) 
an, jedoch wird y; = 1727s; d.h. ys =—1 benutzt und in Heaviside-Ein- 
heiten mit den Kopplungskon- 

g 
stanten B= 4m 
rechnet. Als Energieeinheit ist 
in den Endergebnissen die 
Nukleonmasse M = 1 gesetzt. 


II. Bestimmung derÜbergangs- 
wahrscheinlichkeit 


Wir beginnen mit der Be- 
rechnung des S-Matrixelemen- 
tes nach den Methoden von 
Feynman-Dyson. Die mög- 
lichen Graphen sind von 
Brückner?) angegeben wor- 
den. Wir stellen sie in Abb. 1 
nochmals zusammen. Abb. 14 
und B stellen die erste Nähe- 
rung dar, die Matrixelemente 
sind ~ eg; die übrigen Graphen 
geben zu Matrixelementen eg? 
Anlaß und enthalten die zu 
untersuchenden Korrekturen. 

Die Energieimpulsvektoren 
der Anfangsteilchen sind mit 
k, (Photon) und p, (Proton) 

bezeichnet, die der Endteilchen 

Ce & mit q, (Meson) und pj, (Neu- 
~~~~ Meson ron). Für einige invariante 

Abb.1 ----Photon Größen führen wir Abkürzungen 

ein, nämlich A = 2kq, B=2kp, 

C=2kp' = B—A. Die Matrixelenionte lassen sich unmittelbar anschreiben. 
Sie lauten: 


M,=—2ieg&® ysu), ye Ru) 


i £74 sdlSp (2 + 


*) J.C. Ward, Physic. Rev. 84, 897 (1951). 
RP. Feynman, Physic. Rev. 76, 749, 769 (1949). 
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(P+ R—M)* eu 


My, =—2eg ys (B+ My f dl (B—L— My (el) 
[U + 
Mc, = 2egha'y, fdl(el) 
+ ys (P—Q— My 
Mo, =—eg f al M)> e 
(?— Ys U. 

Bei der Berechnung der Korrekturen können wir die Tatsache benutzen, daß 
in unseren Matrixelementen insgesamt keine Ladungsrenormalisierungsterme 
auftreten; diese kompensieren sich gegenseitig*). Es genügt deshalb, die in My, 
und M 5, enthaltenen Masserenormalisierungen zu berücksichtigen. Um durchweg 
konvergente Ausdrücke zu erhalten, verwenden wir wie üblich®)*) Regulatoren, 
so daß ein Integral / als 

m? 
I= lim f I' (m?) d (m?) 
A 
definiert ist. 
d 
. 
(m M oder u; I’ (m?) 
Weiterhin kénnen wir bei allen Integralen den Hauptwert nehmen. Die etwa vor- 
handenen Imaginärteile der Parameterintegrale fallen bei Bildung der Uber- 
gangswahrscheinlichkeit heraus, da zur zweiten Näherung des Wirkungsquer- 
schnitts nur die gemischten Produkte der Matrixelemente erster und zweiter 
Näherung beitragen. Physikalisch bedeutet dies, daß die Dämpfungsglieder (im 
Sinne der Heitlerschen Theorie) keinen Beitrag liefern. Dies gilt allgemein für 
die zweite Näherung zum Wirkungsquerschnitt, wenn das Element der S-Matrix 
für den betreffenden Prozeß symmetrisch ist bezüglich Anfangs- und Endzustand. 

Zunächst zwei Bemerkungen: 

Einmal wird man erwarten, daß die Korrekturen für kleine Energien relativ 
klein sind. Man überzeugt sich nämlich leicht, daß sie für ~ gleich Null im nicht- 
relativistischen Grenzfall (d.h. k, = 0) verschwinden®). Sie werden also an der 
Energieschwelle unseres Erzeugungsprozesses einen Faktor u = 0,15 enthalten. 
Diese Erwartung wird durch die Rechnung bestätigt. 

Zum zweiten stellen wir fest, daß sich die Matrixelemente M 4, und Mz, exakt 
kompensieren®), sofern man die Massenrenormalisierung beachtet. Wir begnügen 
uns damit, dies durch explizite Rechnung nachzuweisen: 

Es ist zu zeigen, daß die folgende Beziehung gilt (d = qg—k) 

(1) 
eu On (9,9) = en SA Sp {ys (Q—M)* ys (2 + (2+ Q— My} 


ig Sp 19, [R— 8 — M— (2+ O— 


5) Dies gilt nicht für alle Prozesse. Z. B. bleiben die Korrekturen zur Meson-Nukleon- 


streuung auch im N. R. mit « = 0 endlich. 
%) Brückner erhielt für diese Matrixelemente einen kleinen aber endlichen Beitrag. 
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Nach Ward?) läßt sich die rechte Seite umformen, so daß (1) übergeht in 


1 
| Go— 40) = G2 +g (1-2). (la) 
0 


Es genügt, (1a) für die Ableitungen 9), = ds ®, zu beweisen. Nun ist 

eu (9; q) 

fa +2qu + 
[2 — M?) [+ 9)? — + 


Mit Anwendung der üblichen Integrationsmethode ergibt sich 


t 
1 At — 
0 0 


Hieraus erhält man zunächst 


1 
(4, q’) = —ı + | dx [M2 — a — x)? 
0 


und weiterhin 


1 
2 
f dz (qo— qa) do (9°, 9°) 


—i(e fae fa (A2— w) (1 — x)? (le) 


[M? +(42— u) 2)? 


1 
7 A t— p? (1 — 
=—tle [oz [dt 
0 0 


Gemäß (la), (1b) und (lc) ist demnach die obige Behauptung verifiziert. 
Wir haben jetzt die vier noch verbliebenen Matrixelemente zu betrachten. 
a) Beiträge der Graphen B, und B,. 


My, + Mz, =—eg y,(B—M) (B—M) e+Jp, (2) 


1 

—2—M)'dl —2-—M 

0 


+p (1—2) 


B, =—+ [de(1—a) (® —M) (1 — 2) — Mz] 
6 


A 
d 4.4 


mit A, = x* + u? (1— x), wobei die Renormalisation der Masse berücksichtigt 
worden ist. 


1 1—z 
dl (el) 
0 0 


J5, (u?) ; (2b) 
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Nach Einsetzen von (2a, b) in (2), Ausführung der Integration nach u? und Ordnen 
erhält man 


Ms, + B, Mp—ieg? B,(w'y, eu) (3) 


mit 


1 
— Ba + (1 —2) 
=f dee in a? + (1 — 2) 


(4) 


1 
2 
(1-2), 2 Bau (1—2) + 
By de x in + u? (1 —2). 
0 
b) Beiträge der Graphen C, und C3. 
Diese Matrixelemente enthalten keine Renormalisierungsterme; sie sind zu- 
sammen endlich. 


Mo,+ Mc,= eg (5) 
1—z 
di(e 1)(¥’—K — L—M)-! dl (el) (XR + &) 
(1—2) 
1 1—z 
dl & —2-— 
1—z 
Fa 


Durch Integration über u? folgt 


Jc, + Je, (5c) 
e A, 2egx(z+ki) t—(1—z)] e [Pr —R (1—t—M] 
und nach einigen Umformungen 


Mo, 
mit 
1 
__1 [des 
= + (1 — 2) 
0 


(7) 


1 
(dx + Cx (1—2) + 
a= 1+ 


Durch Zusammenfassung der ersten Näherung mit den durch (3) und (6) gegebenen 
Korrekturen können wir nun das gesamte Matrixelement folgendermaßen dar- 


\ 
= 
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stellen: 


M = ieguy,{— Ele et (2 + 1) +0,)— 


Zur Bestimmung des Wirkungsquerschnitts benötigen wir den Ausdruck 


= Wi + Ws. | 
Hierbei wird über Spin und Polarisation im Anfangszustand gemittelt sowie über 
den Endspin summiert. In W berücksichtigen wir im Rahmen unserer Näherung 
Glieder W, ~ a ß sowie W, ~ x ß?. Man findet auf die übliche Weise 

geind 


W, = 220 B+ C,—C, +0,—B,) 
[2- 2B, + 
(= (k, q)). 


Die gewonnenen Ergebnisse hängen von vier Parametern B,, B,,C,,C, ab, die 
wir durch Integrale dargestellt haben. Die Integrationen lassen sich geschlossen 
durchführen. Wir verzichten darauf, die etwas umständlichen Ergebnisse anzu- 
schreiben und stellen statt dessen in Abb. 2 den Verlauf dieser Größen in Ab- 
hängigkeit von B bzw. C graphisch dar’). 


Abb. 2. Als Abszisse ist B bzw. C logarithmisch aufgetragen 


7) Für die Matrixelemente an der Energieschwelle, d.h. für B=u(2 + u) folgt 
hieraus 


Ms, + Ma, = 0,58 Ms; Mc, + Me, = 
4x 
Von Brückner (l. c.) wird statt dessen 


Mz, Mz, 0,38 Ms; Me, Me, 115-5 Me 
angegeben. 
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II. Diskussion des Wirkungquerschnitts 


Wir können nunmehr den Einfluß der Korrekturen auf den Wirkungsquer- 
schnitt untersuchen. Im Laboratoriumsystem gilt 


B=2k;C + u: =2q; 2(B—C) = B[C + u2— cosd Y(C + (10) 


sowie für den Wirkungsquerschnitt 


do 1 


dq, W. (11) 
Mittels (10) kann man hieraus den differentiellen Wirkungsquerschnitt = 


bilden. Der totale Wirkungsquerschnitt wäre durch Integration von (11) über die 
bei vorgegebener Photonenergie möglichen Werte der Mesonenergie zu erhalten. 
Während sich o, leicht angeben läßt, ist die Integration für o, außer im extrem 
relativistischen Fall nicht in einfacher Weise durchführbar. Für k,> M erhält 
man 


Demnach haben die Korrekturen im Grenzfall groBer Energien keine Bedeutung 
fiir den totalen Wirkungsquerschnitt. 


7, 0 
92/6 
0 0 
-02ß + + + + + 
016 05 5 50 500 ke 


Abb. 3. Als Abgzisse ist die Photonenergie k, in Einheiten der Nukleonmasse aufgetragen 


In Abb. 3 ist das Verhältnis der differentiellen Wirkungsquerschnitte für d = 
0, z/, und x dargestellt. Die entsprechenden Kurven für kleine Winkel liegen für 
mäßige Energien unterhalb der 0°-Kurve; die Übergangswahrscheinlichkeit W, 
hat ein schmales Maximum bei 0°. Der Einfluß der Korrekturen auf den differen- 
tiellen Wirkungsquerschnitt wächst für # + 0 mit der Energie unbegrenzt an. 
Es ist jedoch zu bedenken, daß bei großen Energien sich praktisch der ganze 
Prozeß in einem sehr engen Winkelbereich um 0° abspielt. Die gleichfalls einge- 
zeichnte Kurve für das Verhältnis der totalen Wirkungsquerschnitte zeigt den 
schon erwähnten Abfall der Korrektur für extrem große Energien. 

Zusammenfassend können wir folgendes feststellen: 

Für Photonenergien bis zu 1000 MeV ist für den differentiellen Wirkungs- 
querschnitt die Korrektur do, < 0,48 - do,, wobei der Einfluß für große Winkel 
am größten ist. 

Für den totalen Wirkungsquerschnitt ist bei allen Energien 0, < 0,1ß-o,; 
im extrem relativistischen Fall ist o, vernachlässigbar klein. 
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Zusatz bei der Korrektur: Inzwischen wurde mir die Untersuchung von Koba, 
Kotani, Nakai (Prog. Theor. Phys. 6, 849 (1951)) bekannt, die sich eben- 
falls mit den Korrekturen zur Photon-Mesonerzeugung beschäftigt. Unsere Er- 
gebnisse (8), (10), 11) stimmen mit den dortigen, Resultaten (6.1) — (6.5) 
und (6.10) — (6.16) überein. Die Zielsetzung der Arbeiten ist jedoch verschieden. — 
Während es hier darauf ankommt, durch Betrachtung der Korrekturen in Ab- 
hängigkeit von Energie und Winkel für eine bestimmte Theorie. einen Einblick 
in die Güte der Störungsrechnung zu gewinnen, wird in der genannten Arbeit 
der Einfluß verschiedener Kopplungen, neutraler Mesonen und der phänomeno- 
logischen Nukleonmomente untersucht, wobei sich die Auswertung auf eine 
Energie beschränkt. 


Jena, Theoret.-Physikal. Institut der Universität. 


(Bei der Redaktion eingegangen am 10. Juni 1952.) 
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Zum Energie-Impulstensor der Maxwell-Minkowskischen 
Elektrodynamik*) 


Von H. Ott 


Inhaltsübersicht 


Die Unsymmetrie der Energiestrom- und Impulskomponenten des Minkowski- 
schen Tensors rührt von einer Unvollständigkeit dieses Tensors her. Vom Stand- 
punkt der Lorentzschen Elektronentheorie ist derselbe durch ..einen additiv 
hinzutretenden Tensor zu ergänzen, der zum Teil ebenfalls noch von den Maxwell- 
schen Feldgrößen abhängt. Durch diese Ergänzung wird gerade die Symmetrie 
der genannten Tensorkomponenten hergestellt; die bisher scheinbar nur für das 
Vakuum gültige fundamentale Beziehung © = [E §] = c? q überträgt sich nun 
auch auf das Maxwellsche Feld ©, $ in der Materie. 


§ 1. Das Problem 


In der Maxwell-Minkowskischen Elektrodynamik ist die richtige Form 
des Energie-Impulstensors immer noch strittig!). Nun hat jüngst M. v. Laue?) 
in der Strahlgeschwindigkeit einer elektromagnetischen Welle ein wichtiges 
Auswahlprinzip fiir die konkurrierenden Tensoren aufgedeckt und sich unter 
diesem Gesichtspunkt in der Neuauflage seines Buchs?) ‚Die Relativitätstheorie‘ 
für den (unsymmetrischen) Tensor von Minkowski entschieden. Die Unsymme- 
trie desselben hat aber ernste Konsequenzen: So würde sich vor allem die von 
Planck als universell geforderte Beziehüng (in praktischen Einheiten) 


S=[EH] =e (1) 


zwischen Energieströmung S und Impulsdichte g vom Vakuum. nicht auf die 
Materie übertragen lassen. Davon würde unter anderem die Berechtigung der 
Photonenvorstellung in optischen Medien unmittelbar betroffen werden. 

Im folgenden soll nun gezeigt werden, daß der Satz von der Strahlgeschwindig- 
keit die universelle Geltung von (1) keineswegs erschüttert, sondern (freilich erst 


in Verbindung mit der Lorentzschen Theorie) geradezu bestätigt. 


§ 2. Impulsdichte und Kriterium von der Strahlgeschwindigkeit 


Dieses Kriterium besagt, daß die Strahlgeschwindigkeit einer elektrischen 
Welle sich einer Lorentztransformation gegenüber verhalten muß wie die Ge- 


*) Auszugsweise vorgetragen bei der Würzburger Tagung der Physik. Gesellschaft 
in Bayern am 14. Juni 1952 

1) A. Sommerfeld, Vorl. über Theoret. Physik, Bd. III, Wiesbaden 1948, S. 293. 

2) M. v. Laue, Z. Physik 128, 389 (1950). 

3) M. v. Laue, Die Relativitätstheorie, Bd. I, Braunschweig 1952. 
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schwindigkeit eines Massenpunktes (Einsteinsches Additionstheorem 
der Geschwindigkeiten); denn „läuft ein Massenpunkt innerhalb eines Strahls 
mit der Strahlgeschwindigkeit mit, so bleibt er dauernd im Licht, und diese Tat- 
sache ist nicht vom Bezugssystem abhängig‘ (Laue?)). Daran ist natürlich nicht 
zu rütteln. Indessen fällt auf, daß man bei der Anwendung dieses Kriteriums auf 
einen Energie-Impulstensor doch eigentlich nur die 4. Zeile des Komponenten- 
schemas benötigt, nämlich nur die Form des Energiestroms & und der Energie- 
dichte W. So gelangt man z. B. in § 18a bei Laue*) ohne Benutzung irgend- 
eines Energie-Impulstensors bis zu den Gl. (18,10) einschließlich. Unter 
diesen Gleichungen besagt nun (18,9), daß ein gewisser, durch (18,8) definierter 
Vektor vA$ dem Einsteinschen Additionstheorem genügt. Und für diesen Vektor 
folgt weiterhin aus (18,10) die Darstellung (in praktischen Einheiten an Stelle 
der Laueschen C.G.S.): 


[ED] = [EIS DJ] = vA(ED)8 (2) 
oder 
[ES] 
(3) 


Bei keinem dieser Beweisschritte mußte man auf die Form irgendeines Energie- 
Impulstensors zurückgreifen; denn die Größen [E H] und (ED) treten in (3) zu- 
nächst rein formal auf, ihre eigentliche physikalische Bedeutung als Energie- 
strömung und Energiedichte bleibt noch offen. 

Den Energie-Impulstensor, und zwar nur seine Werte © und W, benötigt 
man erst dann, wenn man die (durch (19,20) bei Laue?) definierte) Strahlge- 
schwindigkeit 

S 
»= W (4) 


durch die Maxwellschen Feldgrößen ausdrücken will. Mit den Werten S = [E $] 
und W = 4(€D+ HB) des Minkowski-Tensors (s. 32) wird alsdann ww für 
die ebene Welle identisch mit dem Vektor » A$ der Gl. (3) und genügt somit dem 
Einsteinschen Additionsgesetz. 


Nirgends wird im vorigen die Impulsdichte des Tensors gebraucht. Daher 
verwirft das Kriterium von der Strahlgeschwindigkeit wohl jene unter den kon- 
kurrierenden Tensoren, welche in der 4. Zeile vom Minkowskischen abweichen, 
aber es legt die übrigen Komponenten, insbesondere die Impulsdichte und die 
Symmetrie nicht fest. Jeder andere Tensor, dessen 4. Zeile mit dem Minkowski- 
schen übereinstimmt, wäre ebenfalls zulässig, eine Behauptung, die freilich zu- 
nächst als inhaltsleer erscheinen mag: denn es gibt unter den konkurrierenden 


Tensoren der Maxwellschen Theorie keinen, der die 4. Zeile mit dem Minkowski- - 


schen gemein hätte. Die Situation ändert sich aber, wenn man das Problem im 
Lichte der Lorentzschen, statt der Maxwellschen Gleichungen betrachtet. 
Dann tauchen überhaupt lebhafte Zweifel auf, ob der Maxwellsche Formalismus 
imstande ist, das Tensorproblem in befriedigender Weise zu lösen. Jedenfalls 
zeigt sich, daß der Minkowski-Tensor (wie übrigens jeder andere der Maxwell- 
schen Theorie) vom Standpunkt der Lorentzschen Theorie unvollständig 
ist und durch einen Zusatztensor ergänzt werden muß, der zwar im wesentlichen 
vom ,,verborgenen“ Feld der Materie, aber in gewissem Ausmaß doch auch noch 
vom Maxwellschen Feld abhängt. Ob dieser Zusatztensor auf rein phänomenolo- 
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gischem Wege völlig berechnet werden kann, mag dahingestellt bleiben ; aber soviel 
läßt sich zeigen, daß durch ihn der Energiestrom © und die Energiedichte W des 
Minkowski-Tensors nicht beeinflußt, die Impulsdichte aber auf den Wert 


= +6 reduziert wird (Gl. 38 u. 39), womit sich die Beziehung (1) nun vom 


Vakuum auch auf das Maxwellsche Feld in der Materie überträgt. 


$3. „Ausglättung“ der Lorentzschen Gleichungen 


Ein beliebiger isotroper oder anisotroper Körper werde irgendeinem elektro- 
magnetischen Feld ausgesetzt. Um gewisse Schwierigkeiten bei der Anwendung 
des Gaußschen Satzes zu vermeiden, möge die Dichte des Körpers am Rand 
nicht sprunghaft, sondern stetig gegen Null gehen. Wir legen die Lorentz- 
schen Gleichungen (in praktischen Einheiten) zugrunde, wobei die rasch ver- 
änderlichen Lorentzschen Mikrofelder mit kleinen Buchstaben (e, 6 = u, - h usw.) 
bezeichnet werden sollen. In Viererschreibweise hat man mit dem 


Viererkonvektionsstrom s = fo (5) 
und dem Feldtensor f= 16 _ (6) 
die Feldgleichungen: 
fos = Atv (7) 
0=Aiv f*. (8) 


Die Kraftdichte £=oe-+ [ov,b] und Leistungsdichte q = (ov, e) fassen wir 
zur vierdimensionalen Kraftdichte 


i= (9) 
zusammen, die sich als Vektorprodukt aus s und f darstellen läßt: 
1=[sfl. (10) 


Mit (7,8) und der Rechenregel [Ai v ff] = [Ai v f* f*] — Av v [[f f]]*) folgt aus 
(10) der Lorentzsche Energie-Impulstensor: 


Den Lorentzschen Gln. (7, 8) stellen wir die Maxwellschen gegenüber. Mit 


dem Viererstrom S = {j, ic omaxwen} (12) 
dem Feldtensor P= (13) 
dem Verschiebungstensor H={9,—-ic?} (14) 
und dem Magnetisierungstensor M = {M, ic P}, (15) 


4 
4) [LF H]] ist das Tensorprodukt: [[F Hj];x = F,, Hes — (FH), s. Laue*) 
1 
Ss. 79—81. 
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welcher die durch D = e& E+ ¥, B=, (H + M) definierte Polarisation und 
Magnetisierung zusammenfaßt, lauten die Maxwellschen Gl.: 


S=AtvH 
0=AivF*. 


(16) 


(17) 


Hinzu kommt noch die ,,Verkniipfungsgleichung“ 


F (H + M). (18) 


Durch ‚Ausglättung‘“ der Lorentzschen Gl. ergeben sich bekanntlich die Max- 
wellschen®). Diese Ausglättung (künftig durch einen Querstrich gekennzeichnet) 
läßt sich, da sie als Mittelung über ein (kleines) Raum-Zeitvolumen eine lorentz- 
invariante Operation darstellt, unmittelbar an der Viererschreibweise voll- 
ziehen, womit alles an Übersichtlichkeit gewinnt. 


Aus diesem Grunde erlaube man uns, zunächst das wohlbekannte Ausglättungsver- 
fahren der Lorentzschen Gleichungen mit ein paar Zeilen in Viererschreibweise 
einzukleiden: 

Die ausgeglättete Lorentzsche Gl. (8) oder 0 = div f* wird identisch mit (17), 
wenn man setzt 


j=F. (19) 


Dies ist die bekannte Beziehung e = © und b = wy h = &. Alsdann ergibt die ausge- 
glättete Gl. (7) mit Rücksicht auf (19, 18, 16): 


=Aiv(H+ M)=S+JÄivM. (20) 


0 


Dies ist die Beziehung o0 = j + $ + rot M und o 


= Lorentz 


=o —div §. 


“Maxwell 


$ 4. Der makroskopische Energie-Impulstensor 


Wir stellen uns nun auf den Standpunkt, daß die richtige makroskopische 
Kraft- und Leistungsdichte nicht aus den Maxwellschen Gl., sondern aus der 
ausgeglätteten Lorentzschen Kraftdichte (11) berechnet werden muß. 
Das Ergebnis wird unter Umständen vom üblichen Maxwellschen abweichen. 
Zwar gehen, entsprechend dem § 3, die Maxwellschen Gl. durch Ausglättung aus 
den Lorentzschen hervor, aber hierbei handelt es sich um gine Linearform 
der Felder, bei deren Mittelung die mittleren Schwankungen verschwinden. Bei 
den quadratischen Formen, wie der Kraftdichte, ist eine solche Korrespondenz 
zwischen der Maxwellschen Form und der gemittelten Lorentzschen nicht 
ohne weiteres zu erwarten. Beide werden sich um die mittleren Schwankungs- 
quadrate unterscheiden, und es muß uns im folgenden beschäftigen, wieweit 
letztere eine Rolle spielen, bzw. vernachlässigt werden dürfen. 


Wie gesagt, erblicken wir die richtige makroskopische Kraftdichte Z in dem 


gemittelten Ausdruck (11): 
L = l =a — Ai v 5 (21) 


5) R. Becker, Theorie der Elektrizität, Bd. II, Berlin 1933, S. 110. F. Sauter, 
Z. Physik 126, 207 (1949). 
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wofiir wir mit der Abkiirzung 
1 pen 
— = 6 


L=-AivO 


einfach 


schreiben können. Der Tensor © ist der ausgeglättete Energie-Impulstensor des 
Gesamtfelds, einschließlich des ‚verborgenen‘ inneren Feldes der Materie; er 
ist stets symmetrisch. © wäre noch zu ergänzen durch einen ,,mechanischen“ 
(d. h. nicht-elektrodynamischen) Tensor, der zur Stabilität des Körpers notwendig 
ist, den wir aber hier unterdrücken wollen. 

Die in (22) geforderte Mittelung von [[ff]] kann nicht so einfach wie bei den 
linearen Feldgleichungen vollzogen werden, da es sich jetzt um eine Bilinearform 
der Tensorkomponenten handelt, die nicht unabhängig voneinander gemittelt 
werden können. Zur Weiterbehandlung führen wir die Schwankung £ des Feldes 


um seinen makroskopischen Wert f ein: = f —f, oder mit (19): 
f=F+é. (24) 
Setzt man dies in (22) ein, so folgt, da d = 0 wird: 


In dieser Darstellung zerfallt O in einen ,,Maxwellschen* Tensor 7% — = [[FF]], 


der nur von den Maxwellschen Feldern abhängt, und in einen ,,Sehwankungs- 


tensor“ vo — + [ee], in welchem explizit die Maxwellschen Größen nicht 


mehr vorkommen. Implizit hängt aber 7 sehr wohl noch von diesen ab, wie sich 
sogleich zeigen wird und was auch unmittelbar einleuchtet: Denn das makrosko- 
pische Feld, welches dem Körper durch die Versuchsbedingungen aufgedrückt 
wird, überlagert sich ja nicht einfach additiv dem von Haus aus in der Materie 
vorhandenen Mikrofeld &, sondern modifiziert dieses. Darauf beruht doch gerade 
die Erscheinung der Polarisation und Magnetisierung. 

Daß in der Tat der Tensor T° für sich allein (ohne 7°) nicht als ER 
des Maxwellschen Feldes gelten kann (und ebensowenig der in ihm enthaltene 
dreidimensionale Spannungstensor, auf den man freilich manchmal in elementaren 
Darstellungen der Maxwellschen Spannungen stößt®), vgl. dazu Gl. (26)), lehrt 
ein Blick auf die 4. Zeile des Komponentenschemas 

{Ti T33} = = — -[e =|: Tu, = —- W= - + 
welches nicht nur dem Kriterium von der Strahlgeschwindigkeit (vgl. 3 u. 4), 
sondern unter Umständen sogar dem Energieprinzip widerspricht. Letzteres 
ergibt sich sofort aus der Diskussion der aus 7 resultierenden Leistungsdichte, der 
wir uns nun zuwenden wollen und die uns zugleich einen Hinweis geben wird, wie 
der Gesamttensor © zweckmäßiger in einen ,,Makro‘- und ,,Mikrotensor zu 
zerlegen wäre, 


6) Übersicht über die konkurrierenden Tensoren bei E M ER Die mathem. 
Hilfsmittel des en Berlin 1936, S. 259 u. 277. 
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§ 5. Sachgemäßere Aufspaltung des Gesamttensors 0 


Um die Kraftdichte (23) durch die Maxwellschen Felder auszudriicken, 
führen wir (25) ein und benutzen die Rechenregel —Ai v [[F F]] = [Ai vF, F} 
(da Atv F* = 0). Mit Rücksicht auf (18,16) läßt sich dann weiter schreiben: 


[4iv ]—dive — [Aiv(H + M),F]- Aivw 
= [SF] + [div M, F]—Aivr. 


Unter Benutzung der Werte (12, 15, 13) folgt daraus für die räumliche Kraft- 
dichte L,: 


L, = —Aiv,O = {omaxsw.€ + [7B] — div E+ [B+ rot M, B]}—Atv,r° (26) 
und fiir die Zeitkomponente L, = i/e x Leistungsdichte: 


Ly = — = + (B+ rot M, ©} — Ai m. (27) 


Hierin ist der Beitrag des Tensors 7° in geschweifte Klammern zusammengefaßt. 
Sein Beitrag zur räumlichen Kraftdichte (26), in welchem die Kraft auf die wahre 
Ladung und den Leiterstrom noch ergänzt ist durch eine an der freien Ladung und 
am Polarisations-Magnetisierungsstrom angreifende Kraft, mag vielleicht nicht 
unmittelbar zur Kritik herausfordern. Nichts darf uns aber darüber hinweg- 
täuschen, daß diese Form der Maxwellschen Kraftdichte, die man freilich zu- 
weilen auch in elementaren Darstellungen antrifft®), ohne das Zusatzglied —Ai v, 1? 
unzulässig ist. Dies folgt nämlich sofort aus der Leistungsdichte (27), in der 
neben der Leistung des Leiterstroms j nun konsequenterweise auch eine solche 
des Polarisations- und Magnetisierungsstroms auftritt, und diese Form der Lei- 
stung ohne das Glied — Ai x, r° widerspricht gerade bei den einfachsten dielektri- 
schen Stoffen der Erfahrung. Denn spezialisieren wir auf Stoffe mit temperatur- 
unabhängiger DK und Permeabilitét ohne dielektrische und magnetische 
Verluste”), so hat man doch erfahrungsgemäß als einzig wirkliche Leistung 
des Feldes nur die des Leiterstroms j, welche dem Feld in Form von Wärme 
verloren geht, dagegen keine entsprechende Leistung des Polarisations- und Magne- 


tisierungsstroms. Dem aber widerspricht der Term (B + rot M, €) in (27), der 
im allgemeinen von Null verschieden ist®). Der Widerspruch muß durch den 
Zusatztensor 7° behoben werden, also durch die Bedingung 


(B® + rot M, E) — Ai v, 7° = 0, (28) 


welche dem Umstand Rechnung trägt, daß die Leistung des Polarisations- und 
Magnetisierungsstroms nicht in Wärme umgesetzt, sondern an Ort und Stelle sofort 


7) Letzteres bedeutet: D und B (reelle) lineare Funktionen von € und § mit zeitlich 
konstanten Koeffizienten. 

8) Als einfache Illustration hierzu betrachte man einen vom Gleichstrom j durch- 
flossenen, zylindrischen Draht von der magnetischen Suszeptibilität y. Hier ist zwar 
PB = 0, aber (E rot M) = x (E rot H) = x (j E) ist von Null verschieden und würde eine 
gleichmäßig über den ganzen Querschnitt verteilte Leistung darstellen; zugleich hätte 


man am Drahtmantel eine (negative) Flächenleistung, da M dort einen Flächenwirbel 
besitzt. 
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vom Körper wieder als elektrische Energie aufgenommen wird. Die Trennung 
in die beiden sich kompensierenden Leistungen (28) ist freilich unnatürlich; sie 
rührt von der Zerlegung (25) in die beiden Tensoren 7° und 1° her, die offen- 
sichtlich nicht sachgemäß ist. 

Zu einer zweckdienlicheren Aufspaltung gelangt man aber leicht auf Grund der 
Bedingung (28). Denn diese lautet in Viererschreibweise . 


—Ai-v, [FM] — Aue = 0 (28a) 


und besagt, daß der Tensor 
z=7+[[FM] (29) 


leistungsfrei ist: Adv,r—= 0. Dies (und auch ein Blick auf die Komponente 7",) 
legt die Vermutung nahe, daß in 7° noch ein wesentlicher Anteil des Maxwell- 
schen Feldes steckt, nämlich der „Polarisationsanteil“ [[F M]] mit negativem 
Vorzeichen, den man zweckmäßig heraushebt. Letzteres erreichen wir nun gerade 
durch Übergang zum obigen Tensor r, den wir als besseren Vertreter des ,,ver- 
borgenen“ Mikrofelds ansprechen und demgemäß in die Zerlegung (25) einführen 
wollen. Dann ergibt sich: 


[FM] (30) 


und der hier anfallende Tensor (s. 18) 


7 [tr 0) = ]|— = (31) 


auf den wir ziemlich zwangsläufig geführt wurden, ist nichts anderes, als der 
Minkowskische mit dem leicht verständlichen Komponentenschema: 


T icg=ic[d 3] 


[69] -W=-ZED+HB 


(32) 


‘Die Aufspaltung (30) in die Tensoren 7’ und r halten wir für sachgemäßer, schon 


aus dem Grund, weil die Gesamtleistung nun durch 7 allein gegeben wird: 


Ly = —Aiv, (7 +2) =— din? = 


$ 6. Die Komponenten von 7 


Nun sind wir rasch am gesteckten Ziel, wenn wir uns weiterhin auf Stoffe mit 
temperaturunabhängigen Beiwerten ¢,, und 4, ohne dielektrische und magnetische 
Verluste beschränken. Aus der Thermodynamik®) wissen wir nämlich, daß im 
Ruhesystem die makroskopische tnd durch —O,, gegebene Energiedichte solcher 
Stoffe im Maxwellschen Feld zerfällt in einen elektromagnetischen Anteil von 
der Form 


U,=4(ED+ HB), (33) 


®») R. Becker, Theorie der Elektrizität, Bd. I, Berlin 1944, S. 229. 
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der also vollständig durch die Maxwellschen Größen erfaßt werden kann, und§ 
in einen von der Temperatur # abhängenden Anteil 
U, = f (34)° 


in welchem die Maxwellschen Felder nicht mehr vorkommen. Es ist also 
= U, + U, =U, + fm). (35), 


Da U, mit der Komponente —7',, übereinstimmt (s. 32), können wir statt (35)% 
auch schreiben: 


Or, Tr f (9). 
Der Vergleich mit (30) oder O = T + zeigt, daß 


Ta = — (9) (36)3 
und somit unabhängig vom Maxwellschen Feld ist und offenbar den von? 
Haus aus in der Substanz vorhandenen, een Anteil der inneren“ 
Energie darstellt, mit anderen Worten: — -# ist der elektromagnetische Anteil 


der Massendichte unseres Stoffes. 

Wenn die Substanz nun von einer ebenen Maxwellschen Welle durcheaallll 
wird, läuft nach dem thermodynamischen Ergebnis (35) nur der vom makrosko-? 
pischen Feld & und § abhängige Anteil U, der inneren Energie (35), nicht aber der 
temperaturabhängige !!) Anteil U,, mit dem Felde mit und gibt, wenn tv die Strahl- 
geschwindigkeit der Welle bedeutet, zum Energiestrom S=-U,w= (ED) mw 
Anlaß. Da mw nach dem Satz von der Strahlgeschwindigkeit identisch mit dem? 3 
Vektor vA $ der Gl. (3) sein muß, hat G notwendig die Form: 


S = (ED) w = (ED) 713 = [EG]. sn 

Diese Energieströmung muß in der Komponente 0,, (« = 1, 2, 3) des Gesamt-" 
tensors enthalten sein, die demnach aufspaltet in h 
Ou = + 42), 

worin gx einen nicht mehr von © und § abhängenden Energiestrom darstellt. t 
Da — [€ 9] aber nach (32) gerade die Komponente 7',, des Minkowski-Tensors > 


ist, kann man für die letzte Gleichung auch schreiben: 
Ou. T ~ Pa 
und der Vergleich mit (30) liefert: Die Komponente 


i 
Ta = < Pa 


10) Streng genommen ist nur ein Teil von U, = f(d) in —O,, enthalten, der andere A 
Teil von U, trifft auf den vernachlässigten ‚‚mechanischen‘‘ Energie-Impulstensor, was 3 
im folgenden aber belanglos ist. 


11) Da voraussetzungsgemäß 3 = 0, bleibt die Temperatur bei 
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konstant, s. $8. 
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jst unabhängig vom Maxwellschen Feld (nicht zu verwechseln mit 74, 
da r unsymmetrisch ist!). Schreiben wird das Ergebnis noch einmal an: 


[E + Tu (= £ x Energiestrom ) (38) 


so gilt wegen der Symmetrie von ©: 


Ou = Hla + Tu (= ie X Impulsdichte) (39) 


mit natürlich ebenfalls feldunabhängigem 1,,. 


Die beiden Gl. (38, 39) enthalten das wichtige Ergebnis: Die durch das 
Maxwellsche Feld erzeugte Energieströmung wird vollständig be- 
sehrieben durch den Poyntingschen Vektor S = [E §], die damit ver- 


bundene Impulsdichte durch den Vektor g = + [69]. Somit wird die 
Beziehung (1) auch fiir das Maxwellsche Feld in der Materie gerettet. 


Die tT («= 1, 2,3) bedeuten die Energieströmung (xi/c), bzw. die Impuls- 
dichte (x 7c), welche das innere Feld der Materie von Haus aus hat. Wahrschein- 
lich sind sie im Ruhesystem in den meisten Fällen gleich Null, was sich aber rein 
phänomenologisch wohl nicht so leicht erhärten läßt. Leicht läßt sich jedoch 
zeigen, daß ihr Raumintegral über den ganzen Rn verschwindet. Da nämlich 


7,, unabhängig vom Maxwellschen Feld und somit — En) > aoe = 0 ist, reduziert sich 


die Gl. (28a) oder —Ai v,t = 0 auf die Aussage, daß der dreidimensionale Vektor 
mit den Komponenten 7,,, T,2, T3 überall quellenfrei ist. Da er außerhalb des 
Körpers verschwindet (dort wird mit & und M auch rt = 0), verschwindet nach 
einem bekannten Satz das Raumintegral über jede Komponente. So muß es auch 
sein, denn das verborgene Feld des Kérpers kann im Ruhesystem keinen Gesamt- 
impuls haben’). 

Nicht alle Komponenten von r sind indessen unabhängig vom Maxwell- 
schen Feld, wie schon ein Blick auf Ta = Ty Pa = Ty — {IE Oh 


c? [D B],} zeigt. Ob sich diese Abhängigkeit für alle t,g auf rein phänomenolo- 
gischem Wege ermitteln läßt, ohne auf die atomistische Struktur der Materie ein- 
zugehen, soll dahingestellt bleiben. Wie dem auch sei, das eine ist jedenfalls klar, 
daß der Minkowski-Tensor 7 für sich allein nicht als Repräsentant 
des Maxwellschen Feldes gelten kann. 


$ 7. Kraftdichte und resultierende Kraft 


Mit dem Ergebnis (39) lauten nun die Komponenten L, (« = 1, 2, 3) der räum- 
lichen Kraftdichte unter der Berücksichtigung, daß t,, unabhängig von der Zeit . 


ist: 
ac 


12) Denkbar wäre jedoch ein geschlossener innerer Energiestrom und ein damit ver- 
bundener Drehimpuls des verborgenen Feldes. 


Ann. Physik. 6. Folge, Bd. 11 3b 


A 
et 
- 
oe 
= 
7 


42 Annalen der Physik. 6. Folge. Band 11. 1952 


oder in Vektorform zusammengefaßt 


a6 
(40) 
Das eine Glied ist die Spannungsreaktion, das andere der Rückstoß bei der Änderung 
des Strahlungsimpulses. 

Um die Kraftdichte durch die Maxwellschen Feldgrößen auszudrücken, 


greift man am besten auf 
L=—Aiv(T+ 1) =—Aiv [[F A] Avr (41) 


zurück. Hierbei muß man jedoch beachten, daß —Ai v [[F H]] nicht ohne weiteres 
in das Vektorprodukt [Ai v H, F] übergeht, denn es gilt für zwei zunächst be- 
liebige Sechservektoren F und H (s. Laue’), $. 80): 


— Ai v, [[F H]] + [Ai F*, H*),— [divH,F], 
Die rechts stehende Summe verschwindet, wenn H eine lineare Funktion von F 
mit konstanten Koeffizienten von besonderer Symmetrie ist, nämlich 


= ys mit = 23. Diese Bedingungen sind für unsere speziellen 
yö 


Tensoren F und H erfüllt, bisaufdie Konstanz der az, die am Rand des Körpers, 
wo die Dichte nach Voraussetzung stetig zu Null abfallen soll, als Ortsfunktionen 
anzusehen sind. Daher verschwindet die Summe in (42) zwar in der Zeitkompo- 
nente k = 4, in allgemeinen aber nicht in den Raumkomponenten k = 1, 2, 3, 

Mit Rücksicht darauf folgt aus (41) und [AivH, F] = [SF] für die räum- 
liche Kraftdichte Z, (bei isotropen Substanzen): 


L, = oMaxw. € + [j B] (Z’grade + H? grad u) —Miv,t (43) 


und fiir die Zeitkomponente wegen Ai v,t = 0: 


L,=—Aiv, = =+ 46. (44) 


Die Zusatzkraft —Ai v,t in (43) läßt sich mit den bisherigen Ergebnissen nicht 
völlig berechnen, da nicht alle t,g genügend bekannt sind. Ein Teilbetrag kann 
jedoch noch durch das Maxwellsche Feld ausgedrückt werden, nämlich (z. B, 
bei der 1-Komponente) der Term 


a 
= ax, (TH + Ta 74) Ox, (Tu T 


Aus den Komponenten (32) ergibt sich dann 
3 3 
— Ai dr, 2 Bl [E Oh ome > a or, 1a (45) 


Bemerkenswert ist, daß durch (45) eine Schwierigkeit beseitigt wird, die in der 
elementaren Theorie der Maxwellschen Spannungen auftaucht. Dort erhält 
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man bekanntlich nur das erste Glied = — 3] des Ausdrucks (45) 2), welches auch 
im Vakuum nicht EEE, und in diesem Fall nach (43) eine unver- 
ständliche ‚bewegende Kraft auf das Vakuum“ von der Größe L, = fo Ho [E H]= 


1 2 
ed 
kürlicher Weise, indem man diese Vakuumkraft von vornherein von der Kraft- 
dichte abzieht (vgl. R. Becker®), 8. 237.) Gerade diese Korrektur stellt sich 
nun durch das zweite Glied von (45) automatisch ein. 


[E H] liefern würde. Dieses Dilemma umgeht man SER in etwas will- 


Die Unbestimmtheit der in (45) noch verbliebenen 7,; stört nicht weiter, wenn 
wir die Resultante über den ganzen Körper K bilden: 


R= [L,dv=— +1) dv. 


Denn ein Teil dieses Raumintegrals läßt sich in ein Flächenintegral über die 
„äußere“, ganz im Vakuum verlaufende Oberfläche F des Körpers verwandeln, 
wo r verschwindet, und da wegen der Symmetrie von © ferner 


é [7 é 
ar, (7 + = dx, + Ta = Ta = 
ist, folgt: 


= E(n)df— J S dv (46) 


Un) bedeutet die aus 7 resultierende, an F angreifende Spannung (n die Flächen- 
normale) und hat die Form 


T(n) = E(D n) + (Bn) (ED + HB), 


worin noch D = ¢, E und B = wy H gesetzt werden darf, da F ganz im Vakuum 
liegt. 


Für ein vollständiges System (d. h. keine Ladungen und Ströme außerhalb des 
Körpers K) kann F ins Unendliche verlegt werden, womit das Fläckenintegral 
verschwindet und das Raumintegral über den ganzen Raum zu erstrecken ist: 


Dies ist der Rückstoß, den der Körper erfährt, wenn sich der Strahlungsimpuls 
in oder außerhalb des Körpers, etwa durch Aussendung eines Lichtsstrahls, ändert 
(Schwerpunktsatz). (46 u. 47) sind identisch mit den (um die Vakuumkraft 
korrigierten!) Formeln der Maxwellschen Theorie. Sie ergeben sich jetzt 
ganz von selbst, ohne daß man erst zuvor mit mehr oder minder guter Be- 
gründung die Vakuumkraft von der Kraftdichte abzuziehen hätte. 


18) Und aus energetischen Betrachtungen noch die Zusatzkraft + grad (m 0) 


(vgl. R. Becker®), S. 92), die wahrscheinlich in der letzten Summe von (45) enthalten ist. 
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§ 8. Verallgemeinerung auf beliebige Stoffe 


Es darf nicht verschwiegen werden, daß einer Verallgemeinerung unseres Be- 
weises auf Stoffe mit beliebigen (temperaturabhängigen) Beiwerten ¢ und u einige, 
wenn auch vielleicht nicht unüberwindliche Schwierigkeiten entgegenstehen. Eine 
derselben liegt darin, daß bei dielektrischen und magnetischen Verlusten oder bei 
elektro- und magnetokalorischen Effekten der Ausdruck (28), also Ai v,r, nicht 
mehr verschwindet, da er in (27) ja nun die zusätzliche Verlustleistung geben muß, 
In der Tat wird z. B. der kalorische Effekt einer parelektrischen Substanz (Suszep- 


tibilität y prop. 1/9) durch den Term (% €) beschrieben, wie man der Thermody- 
namik®) entnimmt, nach der eine Substanz bei isothermer Änderung des Felds 


die Wärme dQ = — 9 % Ba abgibt. Für y = 1/0 wird daraus dQ = yEdE = 
EdP und in der Zeiteinheit > — EP. Wenn nun, wie in diesen und anderen 
Fällen die Leistung des Tensors t nicht verschwindet, scheint die Aufspaltung 
0 = T +1, auf der unser Beweis fußte, nicht mehr sachgemäß zu sein. 

Eine andere Schwierigkeit ergibt sich daraus, daß bei temperaturabhängigen 
Beiwerten e und u die Größe 4 (ED + H B) nicht den elektromagnetischen Bei- 
trag zur inneren Energie (= U,), sondern den zur freien Energie darstellt. 

Diese und andere Gründe vereiteln die schrittweise Übertragung des oben 
gegebenen Beweises auf beliebige Stoffe. 


Würzburg, Physikal. Inst. d. Universität. 


(Bei der Redaktion eirgegangen am 9. Mai 1952.) 
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Zur Struktur der Weißschen Bezirke um Höhlungen 
in Oberflächen von Eisen-Silizium-Kristallen. 


Von Bruno Elschner 


(Mit 4 Abbildungen) 


Vor einigen Jahren untersuchten Williams, Bozorth und Shockley’) 
das Verhalten der Weißschen Bezirke auf Oberflächen ferromagnetischer Sub- 
stanzen. Das von den genannten Forschern angegebene Verfahren gestattet, die 
schon vorher von N&el?) geforderte Struktur der Gebiete einheitlicher Magneti- 
sierung um Hohlräume oder unmagnetische Einschlüsse zu studieren®). 

Bei der Durchführung solcher Versuche gingen wir wie folgt vor: Eine Probe 
aus technischem Siliziumeisen (etwa 3 Gewichtsprozent Silizium) wurde kritisch 
gereckt und mehrere Stunden bei etwa 10C0° C im Wasserstoffstrom geglüht. 


Abb. 1 Abb. 2 
Weißsche Bezirke in der Umgebung von Löchern auf einer Eisen-Silizium-Oberfläche. 
Von A nach B verläuft eine Kristallkorngrenze. Die Aufnahme wurde mit dem Zeifschen 
Metallmikroskop „Neophot“ gemacht. Vergrößerung: 1:160 
1) H. J. Williams, R. M. Bozorth u. W. Shockley, Phys. Rev. 75, 155 (1949). 
2) L. Néel, Cah. de Phys. 25, 21 (1944). 
3) L. F. Bates u. C. D. Mee, Proc. Phys. Soc. (A) 65, 140 (1952). 
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Dabei trat eine starke Kornvergrößerung durch Rekristallisation ein. Die durch- 


schnittliche Korngröße betrug nach dieser Behandlung etwa 5 mm?. Die Ober- 
fläche der Probe wurde nach der von Williams, Bozorth und Shockley an- 
gegebenen Methode zunächst mechanisch und dann elektrolytisch poliert. Man 
kann nun durch Veränderung der Stromdichte erreichen, daß die Oberfläche 
nicht vollkommen glatt poliert wird, sondern daß durch den Poliervorgang einige 
Löcher in die Oberfläche gegraben werden. Abb. 1 zeigt die Struktur der Weiß- 
schen Bezirke in der Umgebung so entstandener Löcher. 


In einer Mitteilung von A. Mager‘) werden für kugelförmige Einschlüsse 
Blochwandstrukturen angegeben, die sich ergeben, wenn man die Bedingung der 


Abb. 3 Abb. 4 


Weißsche Bezirke, die sich nach dem Anbringen des Eindruckes E ausbilden. Ver- 
größerung: 1:450, Diagonale des Eindruckes: 16 - 10°? mm 


Streufeldfreiheit stellt. Es wird dort gefordert, daß sich parabelförmige Wände 
ausbilden. In der Umgebung des Loches L (Abb. 1) kann man nun sehr deutlich 
diese parabelförmigen Wände beobachten. Innerhalb dieser Wände verläuft die 
Magnetisierung zylindersymmetrisch um das Loch herum. In der Arbeit von 
Néel (I. e.) ist eine andere mögliche Anordnung der Magnetisierung um einen kugel- 
förmigen Hohlraum erwähnt. 


Abb. 2 zeigt schematisch die Magnetisierungsrichtung in den einzelnen Be- 
zirken. Dabei ist allerdings nicht zu entscheiden, ob die angegebene Pfeilrichtung 
die richtige ist oder ob man zur Darstellung der res alle 
Pfeile um 180° drehen muB. 


4) A. Mager, Wissenschaftl. Zeitschrift der Friedrich-Schiller-Universität Jena 1, 
85 (1952). 
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Um auch die Bezirksstruktur in der Umgebung von anderen künstlich ange- 
brachten Hohlräumen studieren zu können, wurde mit der Vickers-Diamant- 
Pyramide eines Mikrohärteprüfers nach Hanemann (Hersteller: Firma Zeiß, 
Jena) ein Eindruck in die Oberfläche eines kleinen Kristalles gemacht. Um über- 
sichtliche Muster zu erhalten, wurde dafür ein solches Kristallkorn ausgewählt, 
dessen Würfelebene gegen die Oberfläche der Probe nur ganz schwach geneigt 
war. Abb. 3 zeigt das auf diese Weise erhaltene Muster. In Abb. 4 sind die dazu- 
gehörigen Magnetisierungspfeile eingezeichnet. 

Die Grenzen des Eindruckes erscheinen in Abb. 3 nicht scharf, weil sie durch 
Anhäufung von Fe,0,-Teilchen im Hohlraum verwischt werden. Die durch den 
Eindruck hervorgerufenen plastischen Verformungen und die um die Druck- 
pyramide herum austretende Wulst geben anscheinend Anlaß für kompliziertere 
magnetische Verhältnisse in unmittelbarer Nähe des Hohlraumes. 

Die Arbeit wurde auf Anregung von Herrn Prof. Dr. M. Kersten in Zu- 
sammenarbeit mit dem Forschungsinstitut für magnetische Werkstoffe begonnen 
und wird z. Z. weiter fortgeführt. Der Firma Zeiß (Mikroskop-Labor A) danke 
ich für wertvolle Unterstützung beim Anbringen der Eindrücke. 


Jena, Physikalisches Institut der Universität. 


(Bei der Redaktion eingegangen am 28. April 1952.) 
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Quantentheorie einer neuen Ursache 
von Ferromagnetismus 


(Vorläufige Mitteilung) 
Von G. Heber 


Inhaltsiibersicht 


Es wird berichtet über Rechnurgen, welche die kürzlich vom Verfasser 
publizierten Untersuchungen betr. der Rolle der Wechselwirkurg zwischen 
den 4s- und 3d-Elektronen in Ferrmagnetika in wesentlichen Punkten ergän® 
zen und fortführen. Insbesondere sind die Eigenfunktionen der 3 d-Elektronen 
jetzt auch durch Funktionen vom Hund-Mulliken-Typ approximiert worden 
Die erhalteren Spinwellenenergien wurden zur Bestimmung der spontanen® 
Magnetisierurg und der Gesamterergie auf statistischem Wege benutzt. Es zeigteg 
sich, daß die bekannten 7''?-Gesetze nur noch als Spezialfälle auftreten, eine 
Folgerung, die dem Vergleich mit der Erfahrung zugänglich ist und Hinweiseg 
geben sollte, welcher Wechselwirkungstyp bei einer bestimmten ferromagneti=® 
schen Substanz dominiert. 


In einigen kürzlich erschienenen Arbeiten!) spricht Zener die Vermutung aus/# 
daß der Ferromagnetismus in Fe, Co, Ni und gewissen Legierungen und Ver 
bindungen nicht, wie sonst üblicherweise angenommen, durch positives 3d—3d» : 
Austausch-Integral A hervorgerufen wird, sondern durch Zusammenwirken dere 
3d—4s- und 4s—4s-Wechselwirkung, welche Effekte zusammen eine indirekteym 
spinabhängige 3d—3d-Kopplung ganz entsprechender Art wie die direkte Heid 
senbergsche 3d—3d-Kopplung bewirken. Die Zenerschen diesbezüglichen Ab® 
schätzungen!) führen zu recht befriedigenden Resultaten, müssen aber als “7 
grob bezeichnet werden und bedürfen dringend einer Verbesserung. 

Der Verfasser hat in zwei Veröffentlichungen an anderer Stelle?) gezeigt, daBe 7 
die Ausdehnung der Blochschen Spinwellentheorie auf ein System aus N 34 
und N 4s-Elektronen neben den N Atomriimpfen zu einer charakteristischen Ver-™ 
änderung der Energie der Spinwellen in Abhängigkeit von ihrer Wellenzahl führt 
Diese Veränderung besteht im wesentlichen darin, daß das Spinwellenband nicht 
mehr unmittelbar an die Energie des Zustandes ohne Spinwellen anschließt/# 
sondern eine endliche Liicke zwischen Spinwellenband und diesem Grundzustandyy 
bleibt. Hieraus ließ sich qualitativ das magnetische Verhalten eines Systems er-§ 
schließen, in welchem der neue, indirekte Wechselwirkungsmechanismus die@ 
Hauptursache des.Ferromagnetimus ist. Es zeigte sich, daß der Vergleich mit der 


1) C. Zener, Physic. Rev. 81, 440; 82, 403, 88, 298 (1951); vgl. auch I. Isenberai 
Physik. Rev. 83, 637 (1951). ö 

2) G. Heber, Wissenschaftliche Zeitschr. der Univ. Jena, Math.-Nat.-Reihe, 1. Jahrg. 
79 (1952), Fortsetzung: ebenda im Druck. es 
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Erfahrung fiir einen nur schwachen Beitrag dieses neuen Mechanismus relativ 
gum alten, Heisenbergschen spricht. — Jedoch wurde klargestellt, daß es durch- 
aus auch in Substanzen mit A <0 Ferromagnetismus geben kann, falls nur der 
erwähnte indirekte Kopplungsmechanismus genügend stark ist. 

Die beiden zitierten Veröffentlichungen®) litten vor allem noch darunter, daß 
für die 3d-Elektronen Funktionen vom Heitler-Londonschen’), für die 4s- 
Elektronen aber solche vom Hund-Mullikenschen Typ eingeführt wurden. 
Inzwischen haben wir nach dem Vorbilde von Slater‘) auch Rechnungen ange- 
stellt, in welchen wir die 3d-Elektronen in Hund-Mulliken-Zuständen ein- 
führten. Es zeigte sich, daß dies die wesentlichen Eigenschaften der Spinwellen- 
Energie, wie sie früher berechnet wurden, nicht ändert, wenn man Parameter- 
Verschiebungen zuläßt°). Außerdem konnten wir zeigen, daß auch der Übergang 
zu Hartree- und Hartree-Fock-Funktionen als Ausgangsfunktionen des Ein- 
elektronen-Problems höchstens Parameteränderungen in der Spinwellen-Energie 
hervorruft. 

Demnach sind die wesentlichen Züge der verallgemeinerten Spinwellen-Theorie 
weitgehend unabhängig von dem dort eingenommenen Ausgangspunkt der Stö- 
zungstheorie und man kann jetzt mit größerem Recht als vorher die Statistik 
zur Bestimmung des aus dem gewonnenen Termschema folgenden makroskopischen 
Verhaltens des Systems einsetzen. Wir haben in erster Linie die spontane Magneti- 
sierung M(T) in Abhängigkeit von der Temperatur 7’ bestimmt, aber auch den 
magnetischen Anteil Z,, der Kristall-Energie erhält man auf ganz ähnliche Weise. 
Es zeigt sich, daß das Blochsche 7"/:-Gesetz für die Magnetisierung nur mehr als 
Spezialfall erscheint, während i. A. die Sättigung nicht mehr mit 7”/:, sondern 
mit e®/T erreicht wird. Bei alleiniger indirekter Kopplung ergibt sich ein reines 
¢/T-Gesetz, wie schon früher erwähnt. Entsprechend ergibt sich, daß für den aus 
E,. folgenden magnetischen Anteil c,, der spezifischen Wärme das 7"”:-Gesetz 
nur mehr als spezieller Fall unter einer komplizierteren Gesetzesgruppe enthalten 
ist. — Bemerkenswert erscheint mir, daß der indirekte Kopplungsmechanismus 
auch bei zwei- und eindimensionalen Gittern Ferromagnetismus erzeugt, wobei 
allerdings die Curie-Punkte dieser Gitter bei gleicher Kopplungsstärke gegenüber 
einem dreidimensionalen Gitter stark herabgesetzt sind. 

Bezüglich des Vergleichs mit der Erfahrung können die qualitativen Argu- 
mente der zweiten der zitierten Arbeiten quantitativ bestätigt werden. So ist 
2. B. bei alleiniger indirekter Kopplung schon bei etwa 10% der Curie-Temperatur 
absolute Sättigung zu erwarten, also z. B. bei Ni schon bei etwa 60° K eine ganz 
deutliche Abweichung vom 7"/:-Gesetz. Experimentell aber ist das 7”/:-Gesetz 
bis 20° K hinab gut bestätigt®). Man muß daraus schließen, daß bei Ni und wahr- 
scheinlich auch den anderen einfachen Ferromagnetika der neue Kopplungs- 
mechanismus nur eine zweitrangige Rolle spielt. Doch darf er gewiß nicht 
völlig vernachlässigt werden. Zudem halten wir es für nicht ausgeschlossen, daß 
man noch Substanzen finden wird, bei welchen dieser neue Kopplungstyp eine 
größere und vielleicht gar erstrangige Rolle spielt (Seltene Erden, gewisse Legie- 
Tungen). 


3) Die außerdem nicht exakt orthogonal aufeinander waren. 

4) J. C. Slater, Physic. Rev. 52, 198 (1937). 

5) Die Starke der beiden Kopplungen fiihren wir als Parameter ein, beabsichtigen also 
zunächst nicht, die Austauschintegrale usw. auszurechnen. 

6) A. Fallot, Ann. Physique 6, 305 (1936). 


Ann. Physik. 6. Folge, Bd. 11 
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Es wäre wünschenswert, eine entsprechende Theorie für höhere Temperaturen 
zu entwickeln insbesondere für Temperaturen in der Umgebung des Curie- 
Punktes, um das paramagnetische Verhalten mit zur Bestimmung der beiden 
Parameter heranziehen zu können. — Wollte man allein aus dem Tieftemperatur- 
Verhalten von M(T) das ,,Mischungsverhiltnis‘ der beiden Kopplungstypen 
bestimmen, so müßte man wohl noch genauere Messungen bei tiefen Temperaturen 
anstellen, weil sich das 7"/:-Gesetz von den Nachbargesetzen nur sehr wenig unter- 
scheidet. In dieser Hinsicht wäre eine Unterscheidung wahrscheinlich leichter 
bei dünnen ferromagnetischen Filmen, deren M(T)-Gesetze sich stärker von- 
einander unterscheiden. Doch kommen hier natürlich die großen Schwierigkeiten 
des Experimentierens mit dünnen Filmen hinzu. — Die verallgemeinerte Spin- 
wellen-Theorie ist bisher nur auf monoatomare Schichten, noch nicht auf dickere 
Schichten, ausgedehnt worden. 

Die in der ersten Arbeit!) des Verfassers ausgesprochene Hoffnung betreffs 
eines besseren Verständnisses der Größe der Sättigungsmagnetisierung ist leider 
deshalb nicht in Erfüllung gegangen, weil der Anteil der indirekten Wechsel- 
wirkung in den zur Diskussion stehenden Ferromagnetika sehr klein ist und deshalb 
auch die 4s-Elektronen praktisch nichts zu dieser Magnetisierung beitragen können, 
Man wird doch wieder darauf angewiesen sein, entweder eine Kontinuums-Theorie 
der Spindichte zu entwickeln, oder zu versuchen, mit mehreren teilweise besetzten 
3d-Bändern zu arbeiten, um der beobachteten spontanen Magnetisierung der 
Ferromagnetika Rechnung zu tragen. 

Der Verfasser ist bemüht, die Theorie in den angedeuteten Richtungen weiter 
zu entwickeln. 


Jena, Theoretisch-Physikalisches Institut der Universität. 


(Bei der Redaktion eingegangen am 28. Mai 1952.) 
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Zur Theorie der Leitfähigkeit starker nicht assoziierender 
Elektrolyte bei höheren Konzentrationen 


Von H. Falkenhagen, M. Leist und G. Kelbg 
(Mit 3 Abbildungen) 


Inhaltsübersicht 


Im folgenden ist die Theorie der Leitfähigkeit und der Dielektrizitätskonstan- 
ten starker 1—1-wertiger Elektrolyte weiterentwickelt worden in Anlehnung an 
die Theorie von Onsager!), Debye und Falkenhagen?). Diese Weiterent- 
wicklung geschieht unter Benutzung einer kürzlich von Eigen und Wicke®) 
aufgestellten neuen Verteilungsfunktion, welche die Hydratation der Ionen mit 
einschließt, und unter Berücksichtigung des endlichen Ionendurchmessers, d.h. 
eines mittleren Abstandes, bis auf den sich die Mittelpunkte zweier Ionen nähern 


können. Diese neue Verteilungsfunktion nid, welche die Wahrscheinlichkeit in 
der Umgebung eines :-Ions an der Stelle P ein j-Ion an der Stelle Q aufzufinden 
angibt, enthält als neue charakteristische Größe die Besetzungsdichte N, der j- 
Ionen im Sinne von Eigen und Wicke. Letztere haben ihre Theorie auf die 
Thermodynamik der Aktivitäten angewandt; das neue Gesetz für die Aktivitäts- 
koeffizienten in Abhängigkeit von der Konzentration wurde von ihnen an den 
1—l-wertigen Elektrolyten NaCl, KCl, RbCl, CsCl geprüft und in quantitativer 
Übereinstimmung mit den experimentellen Daten gefunden. Die sinngemäße 
Erweiterung dieser Theorie, sowie ihre Anwendung auf die Leitfähigkeit der 
Elektrolyte NaCl, KBr sowie RbCl im stationären Fall und ihre Prüfung an 
Hand der Meßergebnisse ist in dieser Arbeit durchgeführt. Für diese speziellen 
Fälle stimmt die Theorie etwa bis zu der Konzentration von 1 Mol pro Liter Lösung 
mit dem experimentellen Befund überein. 


1. Die allgemeinen Grundgleichungen für nichtstationäre Vorgänge 


In einem starken I eoeits seien Ionen der Sorten 1 ...7...8 mit den je- 
weiligen Ladungen e,...e,...e, vorhanden. Im cm? sind dann hye Zeitmittel 


entsprechend n,...”;...Ms pans verteilt. Bedeuten ni$% und nid Verteilungs- 
funktionen, so gilt fiir nichtstationäre folgende Kontinuitätsgleichung®) 


divg (wi? vi?) divp (wig vig) +2 ot 2 wig) = = 0, (1) 
1) H. Falkenhagen, Elektrolyte (1932) S. 165ff; Electrolytes, International Series 
of Monographs, S. 165ff. 

2) P. Debye u. H. Falkenhagen, Physik. Z. 25, 401 (1928). 

3) M. Eigen u. E. Wicke, Die Naturwissenschaften 38, 453 (1951). Siehe auch 
M. Dutta u. S. N. Bagchi, Ind. J. Phys. 24, 2 (1950); S. N. Bagchi, J. Ind. Chem. 
Soc. 27, 4 (1950). 

7" Näheres siehe unter 1) in der zitierten Monographie Gl. (389). 
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wobei die Abkiirzungen benutzt werden 


n; ni% =w'? und nig =wid. (2) 
vi8 ist die Geschwindigkeit eines j-Ions in der Umgebung eines 7-Ions, ent- 
sprechend ist vid die Geschwindigkeit eines i-Ions in der Umgebung eines j-Ions. 
Zur Anwendung von Gl. (1) auf die Leitfähigkeit müssen wir die Kräfte ein- 
führen, welche die Geschwindigkeit der Ionen verursachen, nämlich elektrische 
Kräfte und Wärmebewegung, da wir von Strömungen in der Flüssigkeit absehen 
können. Die durch Diffusion erzeugte Geschwindigkeit des j-Ions wird jetzt mit 
Hilfe der Verteilungsfunktion von Ei gen und Wicke®°) berechnet. Es ist 


kT j 


| vir) 
ni Ni 


worin o, die Reibungsgröße des 7-Ions bedeutet. Die Verteilungsfunktion wi Q 
genügt der statistischen Formel 
wi? = (4) 


(1-3 7 


wobei N, die Besetzungszahl der j-Ionen ist. Da die Geschwindigkeit des j-Ions 
infolge Einwirkung äußerer Kräfte 


Qj 


ist, erhalten wir insgesamt fiir die Geschwindigkeit 


Q 1 k; 
vip = — (Kip — gradgw 
af wi? ) 
= (5) 
vig = — “Farad puig 
wobei zur Abkiirzung 
k 
ni Nj 
gesetzt ist. 
Die Kontinuitätsgleichung lautet somit 
(wi? —k,T gradgwi}) + (wig — k,T gradp wig) 
1 @ j i 
+ + wig) = 0. (7 


Betrachten wir den Fall eines homogenen stromlosen Elektrolyten, so ergeben 
sich fiir die hier in Frage kommenden Verteilungsfunktionen unter Verwendung 
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0 
e 
der Näherung 


(8) 


Wegen der Symmetrie der Verteilungsfunktionen machen wir fiir die Potentiale 
den Ansatz 


(9) 
P 
ae”) 
und erhalten 
wp = un? {1—( (1- (10) 
Die Funktion g (r) genügt folgender partieller Differentialgleichung 
Ag(r) = g(r) (11) 
mit 
4 N; 
(12) 


Die Lésung dieser Gleichung lautet nach Hiickel5) unter Berücksichtigung eines 
endlichen Ionenradius @ 
“xa -xr 


1 e e 


Somit schreibt sich die Nullösung in der Form 


= (1 r = -\; (wi? = wi). (14) 
Wirkt nun ein elektrisches Feld ©, so sind die Kräfte, welche an den Ionen an- 
greifen 
ai? = — e, gradg + e,€ 


(15) 
=—e, gradpyig +e,€. 


Wir setzen jetzt die allgemeine Verteilung aus der Nullösung und einer Störung 
f,; zusammen. Es ergibt sich dann aus der Kontinuitätsgleichung bei Einführung 
einiger Näherungen und unter Verwendung relativer Koordinaten folgende Dif- 
ferentialgleichung 


(+2 + amy + 2 = Me, (16) 


’) E. Hückel, Physik. Z. 24, 192 (1923). 


<1 folgende Ausdrücke 2 
v;Q N; Vip 
= nn, 1—(1— 
| 
oi P N; 
wg = N,N; 1-(1- = 
iQ 
W; 
(4) 
4 
r)=- 13 
g(r) Dyli+xa r (13) 
| 
d 
(7) = 
7 
eben __ 
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wobei y, = yj + yj angesetzt ist. Ferner bedeutet k,= + . Die Zusatz- 


‘ 
potentiale yj; erfüllen die Poissonschen Differentialgleichungen 


4 
n; Ayj er Ips am 
4m 
n, Ay; ep 


Wir nehmen nun die Zeitabhängigkeit in der Form an f,; = @,,ei®t und erhalten 
unter Beriicksichtigung von (17) fiir die Funktionen G,,; die Differentialgleichungen 

k k 4nn,e 4nn;e . 

T 2) AG, ,— tp Op, 
= 0; ox 


Enthält der Elektrolyt nur zwei Sorten von Ionen, so folgt aus dem System der 
Differentialgleichungen (18) 


= —G,, und = = 0, 


Wir erhalten also für unser Problem nur eine einzige Differentialgleichung 
T =)AG —iwG@ mat) =—(2—“)z @ 
und es gilt die Neutralitätsbedingung n,e, + n,e, = 0. 


Nach Einführung der Relaxationszeit © und der komplexen Größe K lautet 
die Differentialgleichung (19) 


[77 
—— 0 
AG — = — (20) 
& 
Es bedeuten 
(21) 
Q1 Qs 
[7 es ky ky (22) 
1 1 
(23) 
ky ky x 
* 
_ 4% (mei , Me 
Wir führen weiter die Abkürzung ein — 
eı [77 
0: 
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Diese Gleichung geht vermöge der Beziehung 


Ayı = 


in eine Bipotentialgleichung für y; über 


AAy; — K? Ay; = Q (27) 


Nun setzen wir 


D, n 
und es bedeutet 
28) 
EX m 
Die Differentialgleichung für das Zusatzpotential lautet jetzt 
und hat das Integral 
Q —Kr Kr B 
Für unseren Fall muß | 
A,=B,=0und (31) 


sein. Im Inneren der Ionenkugel lautet die Lésung der Potentialgleichung Ayı g=0 
vis = cos Deiet, (32) 


Aus den Randbedingungen an der Oberfläche der Ionenkugel finden wir die Kon- 
stanten B, und C. Es ist 


B, =— (1+ xa)e*4 + (1 + Ka) La: za: (33) 
und 
Q 2 


Das Pctential im Innenraum lautet dann 


. rcos® ett, (34) 


Führen wir eine Funktion x ein, die mit dem frequenzabhängigen Teil der Leit- 
fähigkeit verknüpft ist, so entsteht für das Potential 


rcosdeiet 35 
Yi Xx D,KT cos U € (35) 
und es bedeutet 


= E 
22) 
23 
4) 
= 
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2. Dispersion der Leitfähigkeit und der Dielektrizitätskonstanten 


Ähnlich wie in der Arbeit von Debye-Falkenhagen?) berechnet sich die 
spezifische Leitfähigkeit A zu 


(37) 
mit 
— 


Der Realteil von x ist 


Rey 


zall+za) (41) 
(1-VaR) — ” cos(xa Q) + qwOsin(xa 1-9} 
Es bedeuten R = — 0:6 + @?02+1 und Q= + w?@?— 1. Für 


die finden wir die 


(42) 
mit dem Imaginärteil von x 
mi 
[qwO cos xa VqQ— (1—q) sin xa VaQl— 1409}. 
Durch A = 20) und n = 7 kann auch die molekulare Leitfähigkeit A und 


die molekulare Konzentration y des gelösten Stoffes eingeführt werden. 
Mol 
(Z = 6,02. 10% ) 
3. Leitfähigkeit und Dielektrizitätskonstante bei stationären Feldern 


Um die Leitfähigkeit und die Dielektrizitätskonstante starker Elektrolyte für 
den Spezialfall stationärer Felder zu erhalten, brauchen wir nur in der allgemeinen 
Gleichung w = 0 zu setzen. Es ist 


und 


Damit folgt für die Leitfähigkeit 
A=/o—/,— Au 


(45) 
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mit 
&g 


und für die Dielektrizitätskonstante 


& [0] 
Do -0— D = * 1—qxa(1+xa) 


Die Formeln in sich bei 1—1-wertigen Elektrolyten, wenn wir gleiche 
Besetzungszahlen und gleiche Ionenbeweglichkeiten annehmen. Dann folgt aus 
den Gln. (22), (23) und (24) 


(47) 


| n n 
Für die molare Leitfähigkeit A in praktischen Einheiten finden wir dann 
A = Ag 3D, kT xa(1+xa) N (49 
ne 1000 1 * ) 


Ban 9:10%1l+xa 
und für die Dielektrizitätskonstante 
Do = xa (1-4 V2 xa) )- (50) 


Im Fall verdünnter Lösungen gehen die Gln. (39) und (40) in die von Onsager 
bzw. Debye und Falkenhagen gefundenen Formeln iiber, wie sie in der unter ') 
zitierten Monographie in Gl. (514) bzw. (515) u. (516) und Gl. (539) angegeben sind. 
Denn fiir verdünnte Lösungen können wir den Ionenradius vernachlässigen, d.h. 
a—0 und die Besetzungszahlen werden praktisch unendlich groß, d.h. Noo. 
Der allgemeine Fall von Mischungen starker Elektrolyte läßt sich in Anlehnung 
an Gl. (18) ebenfalls behandeln, ähnlich wie es für den stationären Fall von On- 
sager und Fuoss®) unter einfacheren Voraussetzungen ausgeführt wurde. 


4. Berechnung der Leitfähigkeit der 1-1-wertigen Elektrolyte NaCl, KBr und 
RbCl im stationären Fall 


Bei der Annahme gleicher Besetzungszahlen und gleicher Ionenbeweglich- 
keiten stellen wir die’Leitfähigkeit der Elektrolyte NaCl, KBr und RbCl als Funk- 
tion von der Wurzel aus der Konzentration dar. Wir fassen die Daten für NaCl, 
KBr und RbCl zusammen. 


NaCl: Ag = 108,99 2-1 cm? 
a = 4,2 


na = 1,373 Vy 


60,72 0,203xa__ 
Au = x) 


Ary = 50,50 


6) L. Onsager u. M. Fuoss, J. Phys. Chem. 36, 2689 (1932). 


t 
x 
. 
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KBr: Ago = 132,30 Q71 cm? 


a = 3,8-10-*cm 


xa = 1,243 


81,47 n\2 
Ay, = a_ 1) (1—F) 


Vy n 
An = 5050, 


RbCl: - Ago = 132,6 cm? 
a = 3,2-10-§cm 


Ay, 
xa 
Ay = 48,19 . 


Die Abbildungen 1, 2 und 3 zeigen den Einfluß der Konzentration auf die: 


Leitfähigkeit. Trotz der Annahme gleicher Ionenbeweglichkeiten und gleicher 


iL i iL i i 
0 0b 06 08 10 12 th 
”x x MeBwerte nach Kohlrausch 
---- Theorie unter Berücksichtigung der neuen Verteilungsfunktion 
—— Theorie mit Radiuseffekt mit neuer Verteilungsfunktion 
Onsager -Theorie 


Abb. 1. Molare Leitfähigkeit von NaCl in Abhängigkeit von der Wurzel aus der 
Konzentration bei 18° C 


Besetzungszahlen folgen die Kurven recht gut den Meßwerten von Kohlrausch. 
Wie bei den Aktivitätskoeffizienten in der Arbeit von Eigen und Wicke finden 
wir auch hier eine recht gute Übereinstimmung mit dem Experiment. 
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Die ausführliche Diskussion der in Abschnitt 2 erhaltenen Formeln und die 
Verallgemeinerung für beliebige starke Elektrolyte soll in einer späteren Veröffent- 
lichung vorgenommen werden. Ebenfalls sollen der Wien-Effekt”) und die Vis- 
kosität theoretisch behandelt werden. 
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60+ 
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20} 
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Abb. 2. Molare Leitfähigkeit von KBr in Abhängigkeit von der Wurzel aus der Konzen- 
tration. Bezeichnungen wie in Abb. 1 
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Abb. 3. Molare Leitfähigkeit von RbCl in Abhängigkeit von der Wurzel aus der Konzen- 
tration. Bezeichnungen wie in Abb. 1 
7) Die Berechnung des Wien -Effektes unter Berücksichtigung der endlichen Ionen- 
R größe und der Ionenbesetzungszahlen in Anlehnung an H. Falkenhagen [Physik. Z. 
32, 353 (1931)] ist soeben gelungen. Wir kommen hierauf in Kürze zurück. Auch läßt 
sich die Frequenzabhängigkeit des Wien -Effektes berechnen. Der Einfluß der Ionen- * 
. assoziation wird besonders erörtert. 
ch. 
ion Rostock, Institut für theoretische Physik der Universität. 


(Bei der Redaktion eingegangen am 20. April 1952.) 
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Klassische Statistik 
unter Berücksichtigung des Raumbedarfs der Teilchen 


Von H. Falkenhagen und G. Kelbg 


Inhaltsiibersicht 


M. Eigen und E. Wicke!) haben unter Berücksichtigung des Raumbedarfs 
der Teilchen eine neue Verteilungsfunktion abgeleitet, die bei den Aktivitäts- 
koeffizienten starker Elektrolyte bereits mit Erfolg angewandt wurde. Diese 
Formel gilt jedoch nur für zwei verschiedene Teilchensorten mit gleichen Be- 
setzungszahlen. Im folgenden wird diese Verteilungsfunktion auf beliebig viele 
Teilchensorten unterschiedlicher Größe und somit verschiedener Besetzungs- 
zahlen erweitert. Die Ableitung wird mit Hilfe der Statistik ausgeführt. 


1. Bestimmung der thermodynamischen Wahrscheinlichkeit 


Ein System von @-Teilchen, welche im wesentlichen voneinander unabhängig 


sind, befinde sich in einem Volumen V und habe die Gesamtenergie U. Diese Ge- 


samtheit von G-Teilchen setze sich aus verschiedenen Teilchensorteh. der An- 
zahlen G,,G@,,...G@,, zusammen (@ = 53 a, . Im Raum herrsche ein Kraftfeld, 


welches sich aus einem Potential ¢ herleiten lassen soll. Wir wollen nun die Frage 
untersuchen, in welchem Makrozustand wir ein derartiges System verschiedener 
Teilchen antreffen, d.h. welche Teilchendichten n,(p), nz(9),...n,,(9) sich im 
Gleichgewichtszustand einstellen, wo das Potential y herrscht. Derjenige Makro- 
zustand wird sich ergeben, der die meisten Realisierungsméglichkeiten besitzt. 
Wir haben also das Maximum der thermodynamischen Wahrscheinlichkeit W 
zu bestimmen. Um W auszurechnen, gehen wir so vor, daß wir den ganzen Raum 
in sehr kleine Zellen vom Volumen v,, die wir Elementarzellen nennen wollen, 
unterteilen. Dann durchziehen wir den Raum mit Äquipotentialflächen und fassen 
diejenigen Elementarzellen zu Gruppen zusammen, die zwischen den Äquipotential- 
flächen und y+ Ap liegen. Diese Gruppen werden numeriert. Zur o-ten 
Gruppe mögen z, Elementarzellen gehören, die das Volumen », erfüllen sollen. 
Ferner enthalte diese Gruppe 


Zic Teilchen der Sorte 1 
Z., Teilchen der Sorte 2 


Zmo Teilchen der Sorte m. 


1) M. Eigen u. E. Wicke, Naturwiss. 38, 453 (1951); 89, 108 (1952). Weitere Ar- 
beiten erscheinen in der Z. Elektrochem.; s. a. M. Dutta, Proc. Nat. Inst. Sci. India 18, 
247 (1947); 14, 163 (1948); 17, 27, 445 (1951); Naturwiss. 39, 108 (1952). 
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Da die Teilchen ein nicht. zu vernachlässigendes Volumen besitzen, werden in 
ihnen eine Anzahl dieser kleinen Elementarzellen Platz haben. Diese Platzzahlen 


seien entsprechend 7, 2, - . . P,,, so daß sich also die Teilchenvolumina ausdrücken 
durch 

% = Pi% 

Ve = (1) 


Um = Pm Vo. 
Führen wir noch die Besetzungszahlen N, durch N, = = ein, so erhalten wir die 
i 
Beziehungen 


(i =1,2,...m). (2) 


Es sind in der o-ten Gruppe 


Elementarzellen unbesetzt 


ZısPı Elementarzellen von Teilchen der Sorte 1 besetzt 


ZmoPm . Elementarzellen von Teilchen der Sorte m besetzt. 
Ein Makrozustand in der o-ten Gruppe ist gegeben durch die Zahlen 
Numerieren wir die Elementarzellen durch, so ist ein Mikrozustand dadurch 
festgelegt, daß wir angeben, welche dieser Elementarzellen unbesetzt, welche von 
Teilchen der Sorte 1 usw. besetzt sind. Infolgedessen ergibt sich die Anzahl der 
Realisierungsmöglichkeiten in der o-ten Gruppe zu 


(3) 


Pu Zuo = Ze. (4) 


Multiplizieren wir die Einzelwahrscheinlichkeiten, so erhalten wir fiir W den Aus- 


druck 
m 
| 
W = = (5) 
II Zia! 
0 


W machen wir zu einem Maximum unter folgenden Nebenbedingungen: 


a Exo Zuo =U. (E00 = 9), (6) 


oe 


wobei U die Gesamtenergie ist und daher als konstant angenommen wird. Ein 
Teilchen der Sorte u in der Gruppe a besitzt die potentielle Energie &,o. 
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Die Gesamtteilchenzahl der Sorte u ist konstant: 


= G,. (7) 
Die Anzahl der zur Gruppe o gehörenden Elementarzellen ist konstant: 
PuZuo = Zo. (8) 


2. Aufstellung der Verteilungsfunktion 
Wir bilden den log von W und benutzen die Stirlingsche Formel. Es ist 


floe( 2.2..)' = 1082, (9) 
= b= 
und in der Näherung 
o = = =(0 


Jetzt variieren wir die Z,, und verlangen, daß 6 log W bei den Nebenbedingungen 
verschwindet. Es ergibt sich 


dlogW= > {los ( (11) 
o 4 = 


Nun addieren wir mit geeigneten Multipliktoren ß, x, und A, die Nebenbedingungen 


und setzen die Abänderung 
Das ergibt die Beziehung 
{los ( — los + + Pal =0, (13) 


aus der wir 


m 
log ( 2,2) 10824. + B + + Ac Pu = 0 (14) 
erhalten, und es ist ; 
#3 Z, etn (15) 
Diese Gleichung multiplizieren wir mit p, und summieren über « und bekommen für 
m z, 
= — (16) 


Somit finden wir für die ie 


Z, en 


Zu = (17) 


v=0 
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Beachten wir die Relationen 
1 
N, = 


2 =N,Vo: 


so erhalten wir für die: Teilchendichten 


eP fou 


(18) 


oder 


Nye = ae ‘ (19) 


wobei die Abkiirzung verwendet ist 
Duc = N, (20) 


Die Größen b,. lassen sich aus der Forderung bestimmen, daß für &,, = 0 die 
Verteilungsfunktion in mj, "übergehen muß. Das bedeutet: Es herrscht eine gleich- 
mäßige Verteilung der u-ten Teilchensorte im Raume. 


Aus 
bus 
= Ww=12,...,m) - (21) 
i+ 5 


folgt 
(22) 


Setzen wir noch für ß -- , so erhalten wir für die gesuchte Verteilungs- 
funktion den Ausdruck 


= (23) 
1+ 


Formel (23) geht für m = 2 und gleiche Besetzungszahlen in die von E. Wicke!) 
angegebene über, und die Boltzmannsche Verteilung ergibt sich, wenn die Be- 
0 


N 
setzungszahlen unendlich groß werden, bzw. die Verhältnisse ¥- gegen 1 zu ver- 
v 


nachlässigen sind. 
In der Theorie der Elektrolyte interessiert man sich für die Verteilung der . 


Gegenionen in der Umgebung eines Zentralions. Die Wahrscheinlichkeit wid in 
der Umgebung eines i-Ions an der Stelle P ein j-Ion an der Stelle Q anzutreffen, 
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findet man für Elektrolytmischungen aus Formel (23) zu 


ej ¥ fp 
kr 


iQ e 
1+ 2m. —1 


In dieser Formel bedeuten n, die Teilchenkonzentrationen, N, die Besetzungs- 
dichten, e, die Ladungen der Ionen und yep das von der Ionenwolke der i-Ionen 
am Orte Q erzeugte Potential. Dieses Potential geniigt der Poissonschen Glei- 
chung?) 


kr 
An 2; e 


Ayir=— n; n; Dy m(- . 
+= 


2) Auf weitere Anwendungen dieser allgemeinen Verteilungsfunktion gehen wir 
demnächst näher ein. 


(25) 


Rostock, Institut für theoretische Physik der Universität. 


(Bei der Redaktion eingegangen am 9. Juni 1952.) 
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